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Abstract 


A  typical  lightning  discharge  produces  an  intense  20 gigawatt  electromagnetic  (EM)  pulse 
of  50  to  150  ^s  duration.  A  significant  portion  of  this  energy  propagates  upward  where 
it  interacts  with  the  ambient  electrons  in  the  collisional  lower  ionosphere.  In  the  present 
work  the  Boltzmann  kinetic  equation  for  electrons  and  Maxwell’s  equations  for  the  EM 
fields  are  self-consistently  solved  to  simulate  the  interaction  of  lightning  radiated  EM 
pulses  with  the  lower  ionosphere.  Results  indicate  that  attachment,  ionization,  and  optical 
emissions  produced  in  this  interaction  are  sufficient  to  explain  ‘early’  subionospheric 
VLF  signal  perturbations  and  airglow  brightening  observed  simultaneously  with  lightning. 
This  model  also  provides  the  basis  for  the  quantitative  evaluation  of  a  wide  variety  of 
electrodynamic  and  chemical  processes  in  the  ionosphere  which  occur  in  response  to  the 
energy  released  in  lightning  discharges.  Under  nighttime  conditions,  individual  pulses 
with  peak  amplitudes  of  10-20  V/m  (normalized  to  100  km  free  space  distance)  produce 
changes  in  electron  density  of  1-30%  of  the  ambient  while  a  sequence  of  such  pulses  leads 
to  more  than  100  %  modification  at  altitudes  between  85  and  92  km.  Among  the  optical 
emissions,  the  most  promising  for  observations  are  the  1st  (red)  and  2nd  (blue)  positive 
bands  of  Nj,  which  emit  at  rates  of  7  x  107  and  107  cm-3s-1  at  ~92  km  respectively 
for  a  20  V/m  EM  pulse.  The  height  integrated  intensities  of  the  emissions  reach  4  x  10" 
and  6  x  106  R,  last  for  ~50  //s  and  produce  1,200  and  200  Rayleigh-seconds  for  zenith 
observations. 


Preface 


The  three  years  I  spent  at  Stanford  were  the  most  interesting  and  exciting  years  of  my 
life.  Learning  opportunities  and  means  to  achieve  scientific  results  are  really  superb 
at  this  Farm.  Working  at  the  STAR  Laboratory  fulfilled  my  childhood  and  adolescent 
dreams  about  interstellar  space  exploration.  In  a  sense  I  worked  with  problems  that  are 
not  related  to  stars  themselves  but  I  was  moved  by  the  same  feeling  that  someone  might 
have  when  looking  at  a  mysterious  wonder,  like  a  star  or  the  world  itself. 

However,  my  greatest  admiration  goes  to  the  people  whom  I  met  and  worked  with 
at  Stanford  and  in  the  neighbourhood.  Their  integrity,  generosity,  and  understanding  of 
problems  of  someone  who  arrived  to  a  new  country  were  such  a  cheerful  experience 
during  these  years. 

I  am  very  thankful  to  my  advisor  Prof.  Umran  Inan  and  to  my  associate  advisor  Dr. 
Tim  Bell  for  their  trust  in  me  as  a  potential  student  that  allowed  me  to  cross  the  ocean, 
for  introducing  me  to  a  problem  that  resulted  in  this  thesis,  and  for  their  guidance  in 
solving  it.  I  am  also  thankful  to  them  for  helping  me  solve  numerous  everyday  problems 
which  I  and  my  family  encountered  during  these  three  years.  I  am  very  impressed  by 
their  unique  combination  in  seeing  new  and  exciting  scientific  problems  and  finding  ways 
of  solving  them  in  the  real  world.  Their  enthusiasm  in  going  through  all  the  challenges 
of  scientific  exploration  is  an  example  I  would  like  to  follow. 


PREFACE 


viii 

I  appreciate  collaborations  with  Prof.  Leonid  Kazovsky  from  whom  I  learned  about  the 
importance  of  a  clear  scientific  presentation  and  with  Prof.  Peter  Sturrock  who  showed 
me  the  importance  of  an  ability  for  deep  physical  insight. 

I  am  thankful  to  Professors  Bob  Helliwell,  Martin  Walt,  Ron  Bracewell,  and  Dr.  Vikas 
Sonwalker  for  reviewing  my  papers  before  they  were  submitted  for  publication  and  mak¬ 
ing  a  lot  of  suggestions  for  their  improvements. 

I  would  like  to  express  my  gratitude  to  Professors  Fabian  Pease  and  Greg  Kovacs  for 
taking  time  in  spite  of  their  busy  schedules  and  serving  on  the  committee  for  my  Ph.D. 
defense  and  to  Prof.  Dwight  Nishimura  for  serving  as  the  third  reader  at  a  very  critical 
time  prior  to  the  beginning  of  a  new  academic  year  when  it  seemed  to  me  that  almost 
every  Professor  in  the  Department  was  busy  or  out  of  town. 

I  express  my  delight  in  writing  this  manuscript  using  TeX  macros  written  by  Dr.  Bill 
Burgess  and  I  am  very  thankful  to  him  for  saving  me  a  lot  of  valuable  time. 

I  am  grateful  for  discussions  on  various  topics  and  help  I  got  from  Prof.  Oscar 
Buneman,  Prof.  Don  Carpenter,  Prof.  Tony  Fraser-Smith,  Jerry  Yarbrough,  Bill  Trabucco, 
and  present  and  former  students  of  the  VLF  group:  Lee  Poulsen,  Bill  Burgess,  Sean 
Lev.Tov,  Juan  Rodriguez,  Jasna  Ristic-Djurovic,  Sasha  Draganov,  Viktor  Pasko,  Steve 
Cummer,  Dave  Lauben,  Thomas  Mielke,  Dave  Shafer,  Steve  Reising,  Alan  Brown,  Chris 
Regan,  Mark  Danielson,  Darren  Hakeman,  and  Chem  Lin. 

I  am  grateful  to  Mr.  Tom  Sege  and  Mrs.  Dorothy  Sege,  Mr.  Bob  Thorbum  and  Mrs. 
Jean  Thorbum,  Mr.  Jaroslav  Merkelo  and  Mrs.  Stefanie  Merkelo,  Mr.  Vitaly  Vizyr  and 
Mrs.  Vira  Vizyr,  Mr.  Henry  Milcharski  and  Mrs.  Lyalya  Milcharski,  Mr.  Mike  Sullivan 
and  Mrs.  Oksana  Sullivan  who  generously  supported  me  and  my  family  in  many  different 
ways  during  these  three  years  that  made  my  living  and  working  at  Stanford  much  less 
frustrating  and  more  effective  than  it  would  be  otherwise. 

This  thesis  would  not  be  possible  without  the  experience  and  background  which  I 
have  acquired  working  back  in  Kyiv  and  Moscow  with  Prof.  Nicholai  Kotsarenko,  Dr. 


PREFACE 


ix 


Gennadii  Burlak,  and  Dr.  Vitaly  Chmyrev  whom  I  thank  for  the  time  and  efforts  they 
spent  for  me.  Also  I  am  grateful  to  Dr.  Vitaly  Chmyrev  for  introducing  me  to  Prof. 
Umran  Inan  and  Dr.  Tim  Bell  at  a  crucial  point  at  the  beginning  of  the  long  road  which 
finalized  in  this  Ph.D.  thesis. 

I  want  to  express  many  thanks  to  Mrs.  June  Hua  Wang  and  Mrs.  Jenny  Xu  who  helped 
me  a  lot  with  their  secretarial  work  during  these  years. 

I  dedicate  this  thesis  to  my  mother  Svetlana  and  my  father  Nicholai  Taranenko,  my 
wife  Nina,  my  daughter  Christina,  my  mother-in-law  Raisa  and  my  father-in-law  Leonid 
Zarubin,  and  to  other  members  of  my  family  whose  contribution  is  beyond  words  and 
without  whom  I  would  not  be  what  I  am  now. 

Yuri  N.  Taranenko 

Stanford,  California 
September  25,  1993 


This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under 
grant  NAGW-2871  to  Stanford  University. 


Contents 

Abstract  v 

Preface  vii 

Contents  x 

Tables  xiii 

Illustrations  xiv 

1  Introduction  I 

1 . 1  Characteristics  of  Lightning  Discharges  .  2 

1 .2  Overview  of  Previous  Work  .  4 

1.3  Scientific  Contributions  .  12 

2  Electrodynamics  and  Kinetics  of  the  Interaction  15 

2. 1  Ambient  Profiles  and  Field  Orientations  .  15 

2.2  Description  of  the  EM  Fields  and  the  Electron  Distribution  .  17 

2.3  Time  Evolution  of  the  Electron  Distribution  Function  in  a  Single  Slab  .  21 

2.4  Multiple-slab  Model  of  the  Interaction  as  a  Function  of  Time  and 

Altitude  .  26 

2.5  Representation  of  Injected  EM  Pulses  from  Lightning  .  31 

3  Electron  Density  Changes  33 

3.1  Electron  Density  Changes  for  a ‘Tenuous’ Ambient  Profile  .  34 


x 


CONTENTS 


xi 


3.2  Electron  Density  Changes  for  a  ‘Dense’  Ambient  Profile  .  38 

3.3  Daytime  Conditions  .  39 

3.4  Dependence  on  the  Orientation  of  the  Electric  Field  .  41 

4  Optical  Emissions  43 

4.1  Results  for  a  ‘Tenuous’  Nighttime  D  Region  .  48 

4.2  Results  for  a  ‘Dense’  Nighttime  D  Region  .  53 

4.3  Comparison  with  Earlier  Work  [Taranenko  et  al.,  1992]  .  55 

5  Discussion  59 

5.1  Electron  Density  Changes  and  Associated  Effects  .  59 

5.1.1  Perturbation  of  VLF  Subionospheric  Signals  .  59 

5.1.2  Formation  of  Ionization ‘Bubbles’ over  Thunderstorms  .  62 

5.1.3  Stimulation  of  ULF  waves  .  63 

5.2  Excitation  of  Optical  Emissions  .  64 

5.2.1  Emission  Intensities  Compared  to  Airglow  Background  .  64 

5.2.2  Radiative  Transfer  .  66 

5.2.3  Comparison  with  Boeck  et  al.  [1992]  Observations  .  67 

6  Summary  and  Suggestions  for  Future  Research  7 1 

6.1  Summary  .  71 


CONTENTS 


xii 

6.2  Suggestions  for  Future  Research  .  72 

6.2.1  Dependence  of  the  Lightning  Generated  Electron  Density 

Changes  and  Optical  Emissions  on  the  Temporal  Shape  and 
Duration  of  Lightning  EM  Pulses  .  72 

6.2.2  Application  of  Model  Results  to  Data  on  Subionospheric  VLF 

Signal  Perturbations  .  73 

6.2.3  Three  Dimensional  Modeling  of  the  Interaction  of  Lightning  EM 

Pulses  with  the  Ionosphere  .  74 

6.2.4  Possible  Effects  of  the  Lightning  EM  Pulses  on  the  Thermal  and 

Chemical  Balance  of  the  Upper  Atmosphere  .  75 

6.2.5  Stimulation  of  ULF  waves  .  76 

A  Derivation  of  the  Truncated  Boltzmann  Equations  77 

B  Electron  Impact  Excitation  Cross  Sections  83 

C  Numerical  Algorithms  85 

C.l  The  Boltzmann  Kinetic  Equation  .  85 

C.2  Maxwell  Equations  .  86 

D  Quenching  and  Emission  Rates  89 


Bibliography 


91 


Tables 


D.l  Quenching  and  Emission  Rates 


Illustrations 

1.1  Examples  of  different  electric  field  forms  radiated  by  lightning .  4 

1.2  Examples  of  ‘early’  and  ‘fast’  VLF  signal  perturbations  interpreted  as 

lightning-induced  heating  and  ionization  events  [Inan  et  al.,  1993] .  7 

1.3  Transient  airglow  brighting  observed  from  the  Space  Shuttle  [Boeck  et  al., 

1992] .  9 

2.1  Ambient  profiles  of  electron  and  neutral  densities .  16 

2.2  Electron  loss  function  in  air.  .  20 

2.3  Results  of  a  single  slab  (at  90  km  altitude)  simulation  for  a  fixed  applied 

electric  field  of  E  -  15  V/m  applied  at  time  t  -  0 .  22 

2.4  Results  of  a  single  slab  (at  90  km  altitude)  simulation  for  a  fixed  applied 

electric  field  of  E  =  10  V/m  applied  at  time  t  =  0 .  24 

2.5  Results  of  a  single  slab  (at  90  km  altitude)  simulation  for  a  fixed  applied 

electric  field  of  E  =  5  V/m  applied  at  time  t  =  0 .  25 

2.6  Sketch  of  the  grid  used  to  propagate  EM  pulses  in  the  collisional 

ionosphere .  29 

2.7  Dynamic  test  of  the  numerical  technique .  30 

3.1  Dynamics  of  the  interaction  of  an  EM  pulse  with  the  ionosphere .  35 

3.2  Electron  density  changes .  37 


XIV 


ILLUSTRATIONS  xv 

3.3  Electron  density  changes .  38 

3.4  Electron  density  changes .  39 

3.5  Electron  density  changes .  40 

3.6  Electron  density  changes .  41 

4.1  Energy  levels  and  electronic  band  systems  for  N2  and  Nj .  45 

4.2  Evolution  of  the  interaction  in  space  (altitude)  and  time .  49 

4.3  Intensities  of  optical  emissions .  _5J_ 

4.4  Intensities  of  optical  emissions .  52 

4.5  Maximum  emission  rates  as  a  function  of  the  electric  field  amplitude  jEioo-  53 

4.6  Intensities  of  optical  emissions .  54 

4.7  Maximum  emission  rates  as  a  function  of  the  electric  field  amplitude  jEioo-  56 

5.1  Evolution  of  an  electron  density  disturbance  due  to  successive  EM  pulses 

from  lightning .  63 

B.l  Electron  impact  cross  sections  for  optical  emissions .  84 

B.2  Electron  impact  cross  sections .  84 


1 

Introduction 


On  the  average,  approximately  100  lightning  discharges  occur  per  second  over  our  planet, 
in  approximately  2000  thunderstorm  centers  which  are  active  at  any  given  time  [Cham¬ 
bers,  1967].  The  EM  power  radiated  by  a  cloud-to-ground  return  lightning  stroke  may  on 
the  average  be  20  GW  lasting  for  ~  100  ns  (Section  1.1).  The  cumulative  effects  of  these 
electrical  discharges  influence  many  aspects  of  human  life  and  the  global  atmospheric 
environment.  Among  the  prominent  lightning  effects  which  we  hear  of  in  daily  news 
broadcasts  are  the  disturbances  in  power  and  communication  systems,  ignition  of  forest 
fires,  and  even  sometimes  direct  injuries  and  death  caused  by  cloud-to-ground  lightning 
discharges.  Lightning  also  influences  atmospheric  electric  currents  [e.g.,  Tzur  and  Roble, 
1985;  Uman,  1987;  p.  8],  charged  particle  populations  in  the  magnetosphere  [e.g.,  Voss 
et  al.,  1984;  Burgess  and  Inan.,  1993],  and  the  radio  environment  in  the  Earth-ionosphere 
waveguide  [e.g.,  Uman,  1987;  p.  110;  Lewinstein.,  1976]. 

The  purpose  of  this  research  is  to  consider  a  new  form  of  electrodynamic,  thermal  and 
aeronomic  coupling  of  electromagnetic  (EM)  energy  originating  in  lightning  discharges  to 
the  mesosphere  and  the  lower  ionosphere.  This  new  means  of  coupling  was  discovered  in 
the  course  of  recent  measurements  of  ionospheric  disturbances  associated  with  lightning- 
induced  electron  precipitation  (LEP)  from  the  magnetosphere  [Inan  et  al.,  1988].  The 
latter  is  now’  a  well  known  means  by  which  the  lower  ionosphere  is  regularly  disturbed 
(secondary  ionization,  x-rays,  light  emissions,  and  heating)  as  energetic  radiation  belt 
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electrons  precipitate  as  a  result  of  magnetospheric  interactions  with  lightning-generated 
whistler  waves  [ Irian  et  al.,  1990,  and  references  therein].  The  new  form  of  direct 
coupling  considered  in  the  present  work  is  a  completely  different  process  and  involves 
the  intense  heating  of  the  lower  ionospheric  electrons  by  EM  fields  produced  by  lightning 
discharges  which  in  turn  leads  to  the  production  of  secondary  ionization  and  optical 
emissions  [Irian  et  al.,  1991;  Taranenko  et  al.,  1992]. 


1.1  CHARACTERISTICS  OF  LIGHTNING  DISCHARGES 

In  this  section  we  review  the  known  parameters  of  EM  pulses  from  lightning  discharges 
that  are  important  in  the  course  of  our  study.  More  details  on  general  characteristics  and 
phenomenology  of  lightning  discharges  can  be  found  in  Uman  [1987]. 

Lightning  is  a  transient  electrical  discharge  which  extends  through  the  atmosphere  for 
several  kilometers.  The  most  common  cause  of  lightning  is  the  separation  of  electrical 
charges  associated  with  rainfall.  There  are  different  types  of  lightning,  e.g.,  discharges 
within  the  thunderstorm  cloud  (intracloud  (IC)  discharges),  commonly  seen  cloud-to- 
ground  (CG)  discharges,  and  cloud-to-cloud  (CC)  lightning.  Obviously,  the  directional 
pattern  of  the  radiated  EM  fields  depends  on  the  orientation  of  the  discharge.  In  our  model 
we  consider  the  one  dimensional  (upward)  propagation  of  EM  pulses  and  restrict  ourselves 
only  to  general  comments  on  the  dependence  of  our  results  on  the  orientation  and  spatial 
structure  of  the  discharge.  It  is  important  to  note  that  our  one  dimensional  model  does 
nevertheless  quantitatively  account  for  all  electrodynamic  aspects  of  the  interaction  of 
the  lightning  EM  pulses  with  the  ionosphere.  As  discussed  later  in  Section  5.1.1,  3-D 
modeling  would  be  needed  to  determine  the  horizontal  extent  of  the  disturbed  ionospheric 
region  as  well  as  its  location  with  respect  to  the  causative  IC  or  CG  discharges. 

A  typical  CG  lightning  flash  that  lowers  negative  electric  charge  to  the  ground  has 
a  time  duration  of  about  half  a  second  and  consists  of  several  stages  [Uman,  1987;  p. 
10]:  the  preliminary  breakdown  occurs  within  the  cloud,  the  stepped  leader  initiates  the 


INTRODUCTION 


3 


first  return  stroke  in  a  flash  by  propagating  from  cloud  to  ground  in  a  series  of  discrete 
steps;  the  attachment  process  of  the  stepped  leader  begins  when  the  negatively  charged 
leader  tip  nears  the  ground  and  the  electric  field  terminating  on  sharp  objects  exceeds  the 
breakdown  value  of  air  and  one  or  more  upward-moving  discharges  are  initiated  from 
the  objects;  after  the  attachment  process  the  first  return  stroke  occurs,  which  is  a  ground 
potential  wave  that  propagates  continuously  upward  along  the  previously  ionized  and 
charged  leader  path  (channel),  producing  a  peak  current  near  the  ground  of  typically  30 
kA.  Successive  return  strokes  of  lower  intensities  may  occur  if,  after  the  first  return  stroke, 
additional  charge  is  available  to  the  top  of  the  channel.  The  time  between  successive 
return  strokes  in  a  flash  is  usually  several  tens  of  milliseconds. 

Although  our  model  of  coupling  of  EM  pulses  from  lightning  to  the  lower  ionosphere 
is  applicable  to  different  types  of  lightning  discharges,  the  pulse  shapes  and  intensities 
used  in  this  work  are  mostly  representative  of  EM  pulses  originating  during  a  lightning 
return  stroke.  These  EM  fields  have  an  oscillatory  character  and  typical  durations  of 
50-150  (is  with  peak  electric  field  amplitudes  of  up  to  50  V/m  [Uman,  1987;  p.  110] 
measured  at  horizontal  distances  of  100  km  from  the  flash  (designated  here  as  E'ioo)- 
In  two  series  of  comprehensive  measurements,  the  mean  amplitude  for  CG  flashes  was 
found  to  be  .Eioo  ~  8.8  to  11.2  V/m,  with  amplitudes  greater  than  20  V/m  occuring  ~3  to 
8  %  of  the  time  [Krider  and  Guo ,  1983],  Examples  of  electric  field  waveforms  radiated 
by  lightning  return  strokes  are  shown  in  Figure  1.1. 

The  less  than  150  (is  duration  for  the  first  return  stroke  is  short  in  comparison  with  the 
total  duration  of  a  complete  lightning  flash,  which  is  closer  to  ~500  ms  [Uman,  1987; 
p.  13].  However,  the  largest  electric  current,  EM  power,  and  brightest  optical  emissions 
are  generated  by  the  lightning  channel  during  the  time  period  of  the  first  return  stroke 
[Uman,  1987;  p.  110]. 
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Fig.  1.1.  Examples  of  different  electric  field  forms  radiated  by  lightning.  The  upper 
panel  shows  the  electric  field  of  a  positive  CG  discharge  [Brook  et  al.,  1989],  the  bottom 
left  panel  shows  the  electric  field  of  a  negative  CG  discharge  [Lin  et  al.,  1979]  (the 
scale  is  50  /is  per  division,  the  solid  line  corresponds  to  the  first  and  dotted  line  to  the 
subsequent  strokes),  the  bottom  right  panel  shows  the  electric  field  of  an  IC  discharge 
[  We  id  man  and  K rider,  1979]. 


1.2  OVERVIEW  OF  PREVIOUS  WORK 

If  an  electric  field  of  E\qq  >16  V/m  amplitude  were  to  reach  the  ionosphere  at  100 
km  altitude,  where  the  energetic  electron  mean  free  path  is  about  a  meter  [e.g.,  Rees , 
1989;  Appendix  1,  and  p.  273]  and  the  ionization  potential  is  about  16  eV  [e.g.,  Rees, 
1989;  p.  273],  it  would  cause  avalanche  ionization  of  neutrals.  This  simple  consideration 
illustrates  the  potential  for  direct  ionization  of  the  lower  ionosphere  and  stimulation  of 
optical  emissions  by  lightning  radiated  fields.  On  the  other  hand,  it  is  important  to  note 
that  the  EM  pulse  would  be  attenuated  on  its  way  upward  to  100  km  altitude  as  a  result 
of  its  interaction  with  the  increasingly  dense  electron  population  so  that  the  quantitative 
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consequences  of  the  processes  at  the  important  altitudes  of  80  -  100  km  would  depend  on 
the  amount  of  the  radiated  energy  which  reaches  these  altitudes  where  it  can  significantly 
influence  the  electrons.  Thus,  a  self  consistent  evaluation  of  the  electrodynamic  coupling 
of  a  lightning  generated  EM  pulse  to  the  lower  ionosphere,  properly  accounting  for  the 
energy  losses,  the  evolution  of  the  electron  distribution,  and  the  space-time  evolution  of 
the  EM  pulse,  is  needed  in  order  to  make  a  quantitative  assessment  of  the  consequences 
of  the  interaction  of  the  EM  pulse  with  the  lower  ionosphere. 

The  possibility  of  ionization  of  the  lower  ionosphere  by  lightning  was  considered  ear¬ 
lier.  C.  T.  R.  Wilson  1 1925]  evaluated  the  generation  of  electric  discharges  and  ionization 
at  altitudes  of  60  km  stimulated  by  the  electrostatic  dipole  moment  of  thunderstorm 
clouds.  Bailey  and  Martyn  [1934]  made  a  quantitative  evaluation  of  ionization  due  to 
the  electron  acceleration  by  a  strong  EM  wave  in  the  lower  ionosphere,  and  concluded 
that  electron  density  increases  of  up  to  6.5  times  that  of  the  ambient  can  be  produced  by 
EM  fields  from  lightning.  In  their  work,  Bailey  and  Martyn  [1934]  did  not  consider  the 
attenuation  of  the  EM  wave  as  it  propagates  through  the  collisional  and  highly  conductive 
lower  ionosphere  and  showed  that  a  broad  range  of  extra  ionization  would  occur  at  the  E 
region  (- —  1 00  km)  altitudes  due  to  electric  fields  of  0.4  to  15.5  V/m  and  pulse  durations 
up  to  1  ms.  However,  our  results  (Chapter  3)  show  that  the  electric  field  amplitude  at  100 
km  of  a  typical  lightning  EM  pulse  is  not  sufficient  to  produce  ionization  and  that  the 
region  of  perturbed  ionization  lies  at  ~90  km.  We  also  find  that  the  dominant  mechanism 
of  heating  and  ionization  is  the  process  neglected  by  Bailey  and  Martyn  [1934],  namely 
the  absorption  of  the  wave  power  due  to  its  interaction  with  the  collisional  lower  iono¬ 
sphere.  In  spite  of  these  early  theoretical  works  and  some  early  experimental  evidence  (in 
the  form  of  increased  E  region  electron  density  during  thunderstorms)  obtained  at  about 
the  same  time  [e.g.,  Appleton  and  Naismith,  1933],  most  contemporary  researchers  were 
unaware  of  the  possibilities  of  ionization  of  the  lower  ionosphere  produced  by  lightning 
EM  energy  released  in  the  lower  atmosphere.  Techniques  for  electron  kinetic  modeling 
of  the  collisional  ionosphere  were  developed  years  ago  [e.g.,  Gurevich,  1978]  and  have 
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been  extensively  applied  to  VHF  radar  ionization  of  the  upper  atmosphere  and  the  gen¬ 
eration  of  artificial  ionospheric  mirrors  [Short  et  al.,  1990;  Tsang  et  al.,  1991].  However, 
these  techniques  were  not  applied  to  the  interaction  of  EM  pulses  from  lightning  with 
the  lower  ionosphere.  The  application  to  this  problem  of  oversimplified  models  (e.g., 
[Kucherov  and  Nikolaenko ,  1979])  considering  sharp  ionospheric  bounderies  and  a  fluid 
description  of  electrons  resulted  in  underestimations  of  the  resultant  heating  of  electrons 
leading  to  negative  conclusions  about  possible  ionization  by  EM  energy  from  lightning. 

The  first  contemporary  evidence  for  direct  upward  coupling  of  lightning  energy  to 
the  lower  ionosphere  was  in  the  form  of  ‘early’  or  ‘fast’  perturbations  (sometime  called 
‘early’  or  ‘fast’  Trimpi  events)  of  subionospheric  VLF  transmitter  signals.  The  ‘early’ 
events  occur  <  50  ms  after  radio  atmospherics  [Armstrong,  1983]  or  CG  lightning  [Inan 
et  al.,  1988a],  while  the  ‘fast’  events  are  unusually  rapid  (<  50  ms)  amplitude  changes 
[Inan  et  al.,  1988a,  b].  Examples  of  ‘early’/‘fast’  VLF  perturbations,  observed  in  direct 
association  with  lightning,  are  shown  in  Figure  1.2.  Neither  the  ‘early’  nor  the  ‘fast’ 
aspects  of  these  events  are  consistent  with  LEP  events  which  are  known  to  produce 
subionospheric  VLF  perturbations  with  similar  signatures.  VLF  changes  associated  with 
LEP  events  follow  the  causative  lightning  with  a  0.3-1. 5  s  delay  and  onset  durations 
are  typically  0.5-2.0  s  [Inan  et  al.,  1990].  The  lack  of  any  such  measured  delay  in  the 
‘early’  events  is  indicative  of  a  process  in  which  lightning  directly  disturbs  the  ionosphere 
lying  above  [Inan  et  al.,  1991].  The  observed  recovery  times  of  ‘early’/‘fast’  events  are 
10 — 100  s,  indicating  that  a  change  in  the  lower  ionospheric  electron  density  is  probably 
involved.  In  recent  experimental  work  [Inan  et  al.,  1993]  CG  lightning  flashes  correlated 
with  simultaneous  (within  20  ms)  subionospheric  VLF  wave  perturbations  were  found 
to  be  within  <60  km  of  the  perturbed  VLF  Great  Circle  Paths  and  CG  flashes  with 
measured  first  stroke  amplitudes  as  low  as  i?ioo  =  6  V/m  were  associated  with  detectable 
VLF  perturbations. 

The  direct  upward  coupling  of  lightning  energy  to  the  lower  ionosphere  is  also  con¬ 
sistent  with  observations  from  the  Space  Shuttle  [Boeck  et  al.,  1992]  of  transient  airglow 


INTRODUCTION 


7 


•  Transmitter 
■  Receiver 


48.5  kHz  Amplitude  at  Huntsville  25  Jan  90  0200-0245  UT 


Fig.  1.2.  Examples  of  ‘early’  and  ‘fast’  VLF  signal  perturbations  interpreted  as 
lightning-induced  heating  and  ionization  events  [Inan  et  al.,  1993].  The  top  of  the 
figure  shows  the  map  of  the  region  with  transmitter  (48.5  kHz),  receiver  at  Huntsville, 
and  a  thunderstorm  on  the  path  (+’s  indicate  locations  of  individual  CG  flashes).  The 
top  panel  shows  an  example  of  the  observed  VLF  subionospheric  signals.  The  lower  two 
panels  show  details  of  the  signal  perturbations  and  simultaneous  spheric  observations. 
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brightening  over  a  thunderstorm  center,  measured  using  a  television  video  camera  sen¬ 
sitive  to  visible  light  from  360  to  720  nm  with  peak  sensitivity  at  440  nm.  The  authors 
report  that  there  were  no  luminous  discharges  connecting  the  airglow  luminosity  transient 
to  the  thunderstorm  below.  Only  a  single  example  of  airglow  brightening,  reproduced 
in  Figure  1.3,  was  observed  during  the  Mission  STS-41  flight  on  October  7,  1990.  The 
field  of  view  was  scanned  by  the  camera  starting  at  the  bottom  of  each  figure  and  ending 
above  the  airglow  image.  The  frame  rate  resulted  in  a  time  resolution  of  1/60  of  a  second. 

The  first  frame  shows  the  ocean,  the  horizon,  and  the  ambient  airglow  layer  above 
the  horizon.  The  second  frame  shows  a  brightly  glowing  arc  within  the  airglow  region 
above  the  horizon.  According  to  the  description  by  Boeck  et  al.  [1992],  the  region  of 
enhanced  luminosity  was  10  to  20  km  thick  and  about  500  km  in  apparent  horizontal 
extent.  The  airglow  brightened  to  a  level  about  twice  that  of  the  background  airglow. 
At  the  time  of  the  flash,  the  region  where  the  cloud  flash  would  appear  in  the  next  field 
was  already  saturated  by  the  reflection  of  moonlight  from  the  cloud  bank.  However, 
the  authors  identified  the  weak  begining  of  illumination  from  the  lightning  flash  in  this 
frame.  Under  these  conditions,  a  lightning  flash  would  be  recorded  as  an  increase  in 
the  area  that  was  saturated.  Also,  the  time  interval  between  the  recording  of  the  cloud 
and  the  airglow  images  was  about  1  ms,  long  enough  to  allow  a  stroke  to  have  occurred 
before  the  scan  arrived  at  the  airglow  but  after  the  image  of  the  clouds  was  taken.  The 
third  frame  shows  a  very  large  increase  in  the  brightness  of  the  same  lightning  flash 
accompanied  by  a  decrease  in  intensity  of  the  airglow  transient.  The  authors  conclude 
that  the  image  of  the  ‘blooming’  clouds  leaves  the  impression  of  a  discharge  of  large 
horizontal  extent.  The  fourth  frame  shows  the  airglow  transient  diminishing  in  intensity 
while  the  cloud  luminosity  remains  high  but  decreases  from  maximum  intensity.  Boeck 
et  al.  [1992]  concluded  that  the  decay  time  constant  of  the  optical  airglow  transient  was 
of  the  order  of  10  ms,  in  the  range  of  decay  constants  characteristic  of  the  video  camera. 

The  possible  heating  of  ionospheric  electrons  by  VLF  radiation  from  lightning  was 
first  suggested  by  Inan  [1990]  in  a  paper  concerning  a  recent  observation  of  ionospheric 
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Fig.  1.3.  Transient  airglow  blighting  observed  from  the  Space  Shuttle  [Boeck  et  at 
1992].  The  time  sequence  of  frames  is  as  follows:  top  left,  top  right,  bottom  left,  and 
bottom  right  See  discussion  in  text  for  details 
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heating  by  28.5  kHz  waves  from  a  transmitter  with  100  kW  radiated  power.  However,  it 
was  noted  that  an  additional  mechanism  producing  density  enhancements  would  be  needed 
to  explain  the  observed  ‘early’  or  ‘fast’  effects  in  view  of  their  relatively  long  (10-100  s) 
recovery  times  [Inan  1990].  Generalization  of  the  comparatively  simple  model  of  heating 
[Inan,  1990]  to  the  parameters  of  lightning  discharges  and  to  include  the  production  of 
ionization  was  undertaken  by  Inan  et  al.  [1991].  Assuming  a  Maxwellian  distribution  of 
electrons  to  be  maintained  throughout  the  interaction  process,  and  neglecting  the  reflection 
of  the  EM  pulse  from  the  lower  ionosphere  and  processes  such  as  the  attachment  of 
electrons  to  molecular  oxygen,  this  first  study  concluded  that  significant  new  ionization 
can  indeed  be  produced  by  strong  CG  discharges.  According  to  their  results,  lightning- 
induced  electron  heating  in  the  D  region  can  reach  energies  of  several  eV  and  can 
produce  extra  ionization  sufficient  to  cause  the  observed  perturbations  of  subionospheric 
VLF  signals. 

An  additional  effect  of  lightning-induced  heating  of  the  lower  ionsophere  is  the  pro¬ 
duction  of  enhanced  optical  emissions  by  the  heated  electrons.  The  excitation  of  optical 
emissions  should  be  a  natural  consequence  of  lightning  induced  heating  of  ambient  D 
region  electrons  to  energies  of  several  eV  [Inan  et  al.,  1991;  Rodriguez  et  al.,  1992], 
Taranenko  et  al.  [1992]  used  the  temperature  profile  of  heated  electrons  from  [Inan  et 
al.,  1991]  and  assumed  a  Maxwellian  distribution  of  the  heated  electrons  to  consider  the 
intensities  of  optical  emissions  excited  in  the  D  region  by  the  heated  electrons. 

The  results  of  Taranenko  et  al.  [1992]  suggest  that  lightning-induced  heating  of  iono¬ 
spheric  electrons  leads  to  the  excitation  of  emission  lines  from  atomic  oxygen  (O),  molec¬ 
ular  oxygen  (O2),  molecular  nitrogen  (N2),  ions  of  molecular  oxygen  (O2),  and  nitrogen 
(Nj)  with  intensities  ranging  from  10  to  109  R.  For  example,  the  5577  A  line  from  O 
emits  ~60  Rayleigh  (R)  and  lasts  for  ~350  ms  while  the  1st  and  2nd  positive  bands  of 
N2  emit  ~109  R  but  last  only  as  long  as  the  EM  pulse,  namely  ~50-100/is.  The  bright¬ 
est  emissions  from  the  D  region  above  a  thunderstorm  are  expected  to  be  the  allowed 
transitions,  including  the  1st  and  2nd  bands  of  N2  and  the  1st  negative  and  Meinel  bands 
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of  N£.  The  radiation  time  of  these  bands  is  approximately  equal  to  the  duration  of  the 
EM  pulse  (i.e.,  50  to  100  /xs).  For  the  emissions  of  O  and  O2  considered,  the  radiation 
time  is  limited  by  quenching  and  does  not  exceed  a  few  seconds.  Most  of  the  emissions 
appear  to  be  easily  detectable  with  fast  optical  systems  that  have  exposure  times  of  50 
fi s  to  several  seconds  and  use  the  optical  impulse  from  lightning  to  trigger  the  camera. 

Hence,  initial  theoretical  calculations  [Inan  et  ai,  1991;  Taranenko  et  ai,  1992]  assum¬ 
ing  a  Maxwellian  distribution  of  the  heated  electrons  are  in  qualitative  agreement  with  the 
phenomena  of  ‘early’  or  ‘fast’  perturbations  of  subionospheric  VLF  transmitter  signals 
and  the  excitation  of  a  transient  airglow  brightening  in  the  D  region.  The  assumption 
of  a  Maxwellian  distribution  used  in  these  calculations  has  little  effect  on  the  average 
quantities  (e.g.,  energy  of  the  distribution,  conductivity,  etc.).  However,  such  quantities 
as  the  optical-excitation  or  ionization  rates  that  depend  on  the  number  of  electrons  in  the 
high  energy  (>  10  eV)  tail  of  the  distribution  are  strongly  influenced  by  its  exact  form.  A 
Maxwellian  distribution  is  not  typical  for  electrons  in  a  weakly  ionized  plasma  under  the 
influence  of  an  electric  field  [Gurevich,  1978],  and  the  knowledge  of  the  detailed  micro 
physical  processes  of  field-particle  interactions  in  the  collision  dominated  ionospheric 
plasma  is  therefore  essential  for  understanding  the  key  factors  that  determine  the  vertical 
dimension  and  intensity  of  ionization  and  optical  emissions  in  the  heated  region. 

In  this  thesis  ,  we  undertake  a  fully  kinetic  description  of  the  interaction  of  a  lightning 
EM  pulse  with  the  lower  ionosphere,  using  a  self  consistent  solution  of  the  Boltzmann 
equation  for  the  electrons  and  Maxwell’s  equations  for  the  EM  fields. 
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1.3  SCIENTIFIC  CONTRIBUTIONS 

The  purpose  of  this  Ph.D.  dissertation  research  was  to  theoretically  investigate  a  new  form 
of  electrodynamic,  thermal,  and  aeronomic  coupling  of  energy  originating  in  lightning 
discharges  to  the  mesosphere  and  lower  ionosphere.  The  specific  scientific  contributions 
resulting  from  these  efforts  are  presented  in  Chapters  2-4.  A  brief  summary  of  these 
contributions  is  provided  below.  Most  of  the  results  presented  in  this  dissertation  were 
previously  published  in  [Taranenko  et  al.,  1992,  1993a,  6] 

Chapter  2  A  new  model  has  been  developed  which  involves  the  self-consistent  solution 
of  the  Boltzmann  equation  for  electrons  and  Maxwell’s  equations  for  the 
EM  fields  and  which  simulates  the  interaction  with  the  lower  ionosphere  of 
lightning  radiated  50-150  /zs  EM  pulses. 

Chapter  3  A  detailed  quantitative  assessment  was  made  of  electron  density  changes  in 
the  lower  ionosphere  caused  by  heating  of  the  ionospheric  electrons  by  the 
EM  pulses.  Our  results  provide  the  first  detailed  estimates  of  the  absolute 
values  and  altitude  distributions  of  electron  density  changes  for  different 
intensities  of  EM  pulses  and  ambient  electron  density  profiles  and  form  a 
basis  for  quantitative  explanation  of  all  known  features  of  ‘early’  and/or  ‘fast’ 
VLF  perturbation  events  which  occur  nearly  simultaneously  with  lightning. 

Chapter  4  A  detailed  assessment  was  made  of  optical  emissions  excited  by  heating  of 
the  ionospheric  electrons  by  the  EM  pulses.  Our  results  provide  (i)  the  first 
estimates  of  the  absolute  values  of  major  emissions  and  their  altitude  distri¬ 
butions  for  different  intensities  of  EM  pulses  and  ambient  electron  density 
profiles,  and  (ii)  a  quantitative  explanation  of  Space  Shuttle  observations 
[Boeck  et  al.,  1992]  of  an  enhanced  airglow  transient  in  the  D  region.  We 
also  provide  a  detailed  evaluation  of  the  atmospheric  transparency  for  remote 
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sensing  of  the  airglow  transients  for  space-borne,  aircraft-,  and  ground-based 
observations. 


2 

Electrodynamics  and  Kinetics  of  the 
Interaction 


In  this  chapter,  we  describe  the  basic  model  of  the  electrodynamic  interaction  of  an  intense 
EM  pulse  from  lightning  with  the  collisional  electron  gas  in  the  lower  ionosphere. 

2.1  AMBIENT  PROFILES  AND  FIELD  ORIENTATIONS 

The  altitude  range  of  interest  lies  between  70  and  100  km,  for  which  we  use  the  U.S. 
Standard  Atmosphere  [1976]  for  the  altitude  profile  of  ambient  neutrals  and  assume  the 
composition  to  be  ~80%  N2,  ~20  %  O2,  and  a  varying  concentration  of  O  which  is  ~ 
0.0001%  at  80  km  and  ~  5%  at  100  km  in  accordance  with  the  MSIS  neutral  atmosphere 
model  [ Hedin ,  1987].  Our  results  are  only  weakly  dependent  on  the  ambient  neutral 
temperature,  which  is  taken  to  be  250°  K  throughout  the  region.  The  ambient  neutral 
and  electron  density  profiles  are  summarized  in  Figure  2.1.  We  conduct  calculations 
for  the  ambient  electron  density  profiles  taken  from  Barr  et  al.  [1985]  and  represent 
‘tenuous’  nighttime  (a),  ‘intermediate’  or  ‘dense’  nighttime  (b),  and  daytime  (c),  profiles. 
‘Tenuous’  and  ‘dense’  profiles  respectively  represent  somewhat  more  tenuous  and  more 
dense  D  region  than  the  one  used  in  [ Taranenko  et  al.,  1992]  which  is  also  shown  in 
Figure  2.1  for  comparison. 

In  principle,  the  angles  between  the  wave  electric  field  E,  the  ambient  magnetic  field 
D0,  and  the  plane  of  the  ionosphere  can  be  expected  to  influence  the  details  of  the 
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Fig.  2.1.  Ambient  profiles  of  electron  and  neutral  densities.  Electron  density  profile 
(a)  represents  a  ‘tenuous’  nighttime  D  region,  profile  (b)  represents  a  ‘dense’  nighttime 
profile  for  ‘disturbed’  geomagnetic  conditions  and  (c)  represents  daytime  conditions.  The 
electron  density  profile  used  in  our  earlier  paper  [Taranenko  et  at,  1992]  is  also  shown 
for  comparison. 


interaction.  For  simplicity  we  assume  that  the  ambient  magnetic  field  lies  in  the  plane 
of  the  lower  ionospheric  boundary  (i.e.,  horizontal)  and  the  incident  EM  wave  is  planar. 
We  consider  two  extreme  cases  with  E  ||  Bq  and  E  _L  Bq.  With  these  assumptions,  our 
model  results  are  most  directly  applicable  to  the  interaction  of  lightning  EM  fields  with 
the  ionosphere  near  the  geomagnetic  equator.  It  is  shown  in  Section  3.4  by  comparing 
results  for  E  ||  Bq  and  E  L  Bq  that,  due  to  the  intense  nature  of  the  heating  for  typical 
cases,  the  plasma  is  collision  dominated  so  that  the  general  characteristics  of  the  coupling 
do  not  depend  critically  on  the  orientation  of  the  magnetic  field. 
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2.2  DESCRIPTION  OF  THE  EM  FIELDS  AND  THE  ELECTRON 
DISTRIBUTION 


The  propagation  of  EM  fields  is  described  by  the  one  dimensional  (1-D)  Maxwell’s 
equations: 

dBz  1  dEy  _  4tt  dBy  1  dEz  _  4tt 

dx  c  dt  c  y  '  dx  c  dt  c  2  r?  n 

dEz  1  dBy  dEy  1  dBz  A  1  u 

ax  c  dt  dx  +  c  dt  0 

where  c  designates  the  velocity  of  light  in  a  vacuum,  x  is  the  vertical,  2  is  the  horizontal 
axis  directed  along  the  ambient  magnetic  field  B0,  and  y  is  the  horizontal  axis  directed 
in  accordance  with  a  right  handed  x,y,z  coordinate  system.  The  electron  distribution 
function  f(r,v,t )  is  related  to  the  electric  current  J  and  electron  density  Ne  via  the 
subsidiary  equations: 


J 


=-'/ 


vf(r,v,t)dv\  N( 


-J 


f(r ,  v,  t)dv 


(2.2) 


where  e  is  the  electron  charge.  Since,  as  will  be  seen  later  (Section  2.4),  the  random 
electron  velocity  is  much  larger  than  its  average  directional  velocity,  we  expand  /(r,  v,  t ) 
in  spherical  functions  of  zero  order  keeping  only  the  first  two  terms  [Allis,  1956] 


f(r,  v ,  t)  ~  /0(r,  v,  t)  + - + 


(2.3) 


where  fo(r,v,t)  is  the  symmetrical  part,  whereas  lihiLlJl  \s  the  directional  part.  For 
our  case  (2.3)  is  a  good  approximation  as  long  as  e2E2 /mlv2  <C  1  which  is  satisfied 
for  all  parameters  considered  here  (see  Appendix  A),  and  throughout  the  process  of  the 
interaction,  where  m  is  the  electron  mass,  v(e)  is  the  total  effective  electron  collision 
frequency  with  neutrals,  and  e  is  the  average  energy  of  the  electrons.  Substitution  of 
(2.3)  into  (2.2)  gives: 


J  = 


4ne 


f 


Q  — ♦ 

v  f\(r,v,t)dv. 


Ne 


=  4t t  v  fo(r,v,t)dv 

Jo 


3 
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which  explains  the  physical  meaning  of  first  two  expansion  terms  in  (2.3):  fo(r,v,t) 
can  be  used  to  calculate  the  electron  density  whereas  f\(r,  v,t)  is  useful  for  calculating 
the  particle  flow  in  configuration  space.  For  a  weakly  ionized  plasma,  the  Boltzmann 
equation  for  the  electron  distribution  function  can  be  reduced  to  the  following  system  of 
equations  (see  Appendix  A) 


dfo 

dt 


e  d 
3 mv2  dv 


+ 


J_d_ 

2v2  dv 


Tdfo 

m  dv 


(2.4) 


dfl 

dt 


vE  dfp 
m  dv 


Bx/, 


me 


(2.5) 


C  _  qrot  ,  qvib  . 
in  “  °0m  *1"  °0tn  ^ 


sZ * sf" 


,  qatt 
Om  ^  °0 in 


.  Cion 
+  °0*  n 


(2.6) 


The  second  term  on  the  right  hand  side  of  (2.4)  describes  contributions  from  elastic  col¬ 
lisions  of  the  electrons  with  neutrals,  with  6  being  the  fraction  of  electron  energy  lost  per 
collision,  vei(v)  being  the  frequency  of  elastic  collisions,  and  T  being  the  temperature 
of  neutrals.  5b,n  *s  the  part  of  the  Boltzmann  collision  integral  representing  inelastic 
collisions  consisting  of  rotational,  vibrational,  optical,  dissociative,  dissociative  with  at¬ 
tachment,  and  ionization  losses.  The  rotational  term  can  be  represented  as  [Gurevich, 
1978;  p.  53] 


qrot  _  u 

0,n  “  dv  L 


%PNvorot  (Tdfo 


(=$♦•*)] 


(2.7) 


m  \m  dv 

where  (3  is  the  rotational  constant  for  a  specific  molecule  (see  Appendix  B),  N  is  the 
density  of  molecules,  and  <rrot  is  the  rotational  cross  section  (see  Appendix  B).  For  a 
weakly  ionized  plasma  the  vibrational,  optical,  and  dissociative  losses  are  represented  in 
a  similar  manner  [Holstein,  1946] 


Soin  — 


2 N_ 
mv 


^2  [(e  +  e>)/o(e  +  ej)crj(e  +  tj)  -  e/0(e)cr_,(e) 
j 


(2.8) 


2 

where  t  -  is  the  energy  difference  between  the  ground  and  j  levels,  and  cry  is 

the  corresponding  cross  section  (see  Appendix  B  for  numerical  values).  The  attachment 
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term  is  given  as  [Gurevich,  1978;  p.  79] 


Sq"  =  -No2<Tatt(v)vf0 


(2.9) 


where  No2  is  the  density  of  O2  molecules,  and  oatt  is  the  attachment  cross  section  (see 
Appendix  B).  Ionization  processes  are  described  by  the  following  equation  [Gurevich, 
1978;  p.  77] 


(2.10) 


+  —  (*  +  e*)/o(e  +  <0^,o"(e  +  «)  -  e/o(e  V°n(0 


mv 


where  vton  is  the  electron  velocity  corresponding  to  the  ionization  threshold  for  the 
particular  molecule  (O2  or  N2),  vxon(v' ,  v)  is  the  number  of  electron-ion  pairs  produced 


in  one  second  by  an  electron  having  a  velocity  v'  and  leading  to  the  appearance  of  a 


new  electron  having  a  velocity  v,  e*(e)  is  the  average  energy  lost  by  electrons  in  a  single 
ionization  as  given  in  calculations  using  Opal  et  al.  [1971]  data  on  the  secondary  electron 
spectra  and  the  ionization  potential  approximation  used  by  Richards  and  Torr  [1990].  The 
dissociative  attachment  of  electrons  (C>2+e  — »  O  +  O-)  and  ionization  (M2+e  — ►  M2+2e, 
where  M2  designates  a  N2  or  O2  molecule)  are  processes  in  the  model  that  affect  the 
electron  density.  We  neglect  existing  spatial  gradients  of  the  distribution  function  since 
they  are  only  important  in  the  D  region  over  time  intervals  longer  than  several  hundred 
seconds  [Tomko  et  al.,  1980]  and  the  duration  of  the  lightning  EM  pulse  is  typically  <100 
/is.  The  total  electron  energy  loss  function  F,  where 


(2.11) 


(where  the  summation  is  performed  over  all  inelastic  loss  processes  and  C2  is  the  con¬ 
centration  of  the  particular  sort  of  neutrals,  Sej  is  the  energy  loss)  is  shown  as  a  function 
of  energy  in  Figure  2.2  and  is  similar  to  that  given  in  Figure  7  of  Tsang  et  al.  [1991].  We 
note  that  there  are  two  major  energy  barriers.  The  first  one  is  due  to  the  N2  vibrational 
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barrier,  extending  from  1.5  to  3  eV.  Inelastic  losses  in  the  range  3  to  6  eV  are  much 
smaller  than  the  vibrational  losses  of  the  first  barrier.  The  second  barrier  appears  in  the 
range  above  ~6  eV,  and  is  due  to  the  optical,  dissociative,  and  ionizational  excitations 
of  the  ambient  neutrals. 
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Fig.  2.2.  Electron  loss  function  in  air.  Equation  (2.11)  and  numerous  electron  impact 
cross  sections  (Appendix  B)  were  used  to  compile  this  function. 
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2.3  TIME  EVOLUTION  OF  THE  ELECTRON  DISTRIBUTION  FUNCTION 
IN  A  SINGLE  SLAB 

Even  after  the  simplifying  assumptions  introduced  in  the  previous  section,  the  solution  of 
the  one  dimensional  Maxwell’s  equations  (2.1)  simultaneously  with  the  equations  (2.4) 
and  (2.5)  for  the  evolution  of  the  electron  distribution.function  as  given  by  (23)  for 
time  varying  fields  over  the  altitude  range  of  interest  (70  to  100  km)  requires  extensive 
computer  resources.  It  is  thus  useful  to  make  additional  simplifications  by  taking  into 
account  the  parameters  of  our  particular  problem.  For  some  insight,  we  first  examine  the 
time  evolution  of  the  distribution  function  in  a  single  slab  for  a  given  applied  electric 
field. 

We  solve  equations  (2.4)  and  (2.5)  including  all  loss  terms  (2.6)  under  the  assumption  of 
a  constant  electric  field  using  a  numerical  technique  for  the  Bo)t7.mann  transport  equation, 
as  described  for  example  by  Rockwood  and  Greene  [1980],  which  conserves  energy  and 
the  number  of  particles.  Details  of  the  numerical  technique  are  given  in  Appendi^-C. 

The  time  evolution  of  the  distribution  function  in  a  single  slab  located  at  90  km  altitude 
for  a  constant  electric  field  of  15  V/m  applied  for  20  //s  is  shown  in  Figure  2.3  (b). 

The  distribution  function  is  initially  a  Maxwellian  with  temperature  equal  to  that  of 
the  neutrals  (~250°  K).  After  the  onset  of  the  electric  field,  the  energy  Of  the  electrons 
increases,  reaching  a  constant  value  of  ~5.7  eV  in  about  5  ps  (Figure  2.3  (a)),  at  which 
time  the  distribution  function  has  reached  its  quasi-equilibrium  state.  This  state  is  char¬ 
acterized  by  an  almost  flat  distribution  below  ~1.2  eWwith  a  drop  between  1.2  and  3  eV 
due  to  the  vibrational  excitations  of  N2  (see  Figure  2.2).  A  second  flat  region  exists  from 
3  to  6  eV  due  to  the  gap  in  the  electron  loss  function  as  can  be  seen  from  Figure  2.2  and 
as  discussed  above,  followed  by  a  precipitous  drop  as  electrons  reach  energies  above  the 
major  optical,  dissociative,  and  ionization  thresholds  of  the  atmospheric  constituents  and 
the  loss  rate  increases  rapidly  with  energy  as  seen  in  Figure  2.2.  The  corresponding  time 
evolutions  of  the  electric  current  density,  average  energy,  ionization  and  attachment  rates 


Normalized  Rates  (crrf  Js) 
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Fig.  2.3.  Results  of  a  single  slab  (at  90  km  altitude)  simulation  for  a  fixed  applied 
electric  field  of  E  =  15  V/m  applied  at  time  t  =  0.  Top  left:  Time  evolution  of  the 
normalized  electric  current  and  average  electron  energy.  Top  right:  Time  evolution  of 
th^  electron  distribution  function.  The  electron  distribution  function  is*"sfiown  at  0.022 
ps,  0.25  psy  0.75  pSy  2.3  pSy  14  ps>  and  20  ps.  Bottom  left:  Time  evolution  of  the 
attachment  and  ionization  rates.  Bottom  right:  Time  evolution  of  the  optical  emission 
rates  for  selected  lines.  The  current  density  and  excitation  rates  are  normalized  to  one 
electron  per  cubic  centimeter. 
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are  shown  in  Figure  2.3  (a)  and  (c).  It  is  seen  from  Figure  2.3  (c)  that  the  attachment 
rate  reaches  its  equilibrium  value  over  a  shorter  time  interval  than  the  ionization  rate. 
This  behaviour  is  due  to  the  lower  threshold  energy  of  attachment  (~6  eV)  in  compari¬ 
son  with  ionization  (~16  eV).  The  excitation  rates  of  selected  optica]  emissions  shown 
Figure  2.3  (d)  behave  in  a  similar  manner  with  higher  states  (e.g..  the  1st  negative  of  N} 
and  Oj.  and  Meinel  of  N£)  reaching  equilibrium  values  at  slightly  later  times  compared 
to  lines  with  lower  thresholds  (6300  and  5577  .4).  The  rational  for  selection  of  the 
particular  emission  lines  and  bands  is  given  in  Chapter  4. 

Two  other  examples  of  single  slab  calculations  for  applied  electric  fields  of  10  and 
5  V/m  are  shown  in  Figure  2.4  and  2.5.  We  note  that  the  system  reaches  equilibrium 
state  over  somewhat  longer  time  for  smaller  electric  field  amplitudes.  For  example,  90 
%  value  of  the  equilibrium  ionization  rate  (see  Figure  2.3  (c),  2.4  (c),  and  2.5  (c)) 
is  reached  in  ~3  /is  for  E  =  15  V/m,  in  ~4  /is  for  E  =  10  V/m,  and  in  ~15  /is  for 
E  =  5  V/m.  This  behaviour  is  due  to  slower  development  of  the  high  energy  tail  of 
the  distribution  (which  determines  the  ionization  rate)  for  lower  electric  field  amplitudes. 
However,  since  the  electric  current  is  determined  by  the  bulk  of  the  electrons,  the  1 10 
value  (the  current  becomes  smaller  with  time)  of  the  equilibrium  current  is  established 
in  smaller  time  intervals  (Figure  2.3  (a),  2.4  (a),  and  2.5  (a))  ~2  /is  for  E  =  15  V/m. 
~3  /is  for  E  =  10  V/m,  and  ~3  /is  for  E  =  5  V/m.  We  note  the  more  than  three 
orders  of  magnitude  difference  in  the  equilibrium  ionization  rates  for  E  =  15  and  E  =  5 
V/m  which  underscore  the  highly  nonlinear  nature  of  the  interaction  and  the  resulting 
ionization  process. 

As  will  be  seen  from  the  results  presented  in  Chapters  3  and  4,  the  major  electron 
density  changes  and  excitation  of  optical  emissions  for  the  electron  profiles  considered 
in  this  thesis  occur  at  or  below  90  km  altitude.  Hence,  we  have  chosen  to  present  our 
single  slab  simulation  results  for  this  altitude. 


Concerning  other  altitudes  and  electric  field  amplitudes,  we  note  that  from  equation 
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Fig.  2.4.  Results  of  a  single  slab  (at  90 
km  altitude)  simulation  for  a  fixed  ap¬ 
plied  electric  field  of  E  -  10  V/m  ap¬ 
plied  at  time  t  =  0.  Top  left:  Time  evolu¬ 
tion  of  the  normalized  electric  current  and 
average  electron  energy.  Top  right:  Time 
evolution  of  the  electron  distribution  func¬ 
tion.  The  electron  distribution  function  is 
shown  at  0.022  /js,  0.25  ps,  0.75  /is,  2.3 
ps,  14  /is,  and  20  ps.  Bottom  left:  Time 
evolution  of  the  attachment  and  ionization 
rates.  The  current  density  and  excitation 
rates  are  normalized  to  one  electron  per 
cubic  centimeter. 


(C.l)  in  Appendix  C  that  the  first  term  on  the  right  hand  side,  namely  the  driving  force 
Fe  due  to  the  electric  field  E,  is 

Fe~E2/v  (2.12) 

The  diffusion  time  of  the  distribution  function  over  a  unit  of  velocity  or  energy  space 
under  the  influence  of  the  electric  field  is  inversely  proportional  to  this  force  and  is 
estimated  as  ~  v/E1.  The  collision  frequency  v  is  proportional  to  the  neutral  density 
which  changes  by  a  factor  of  ~6  over  a  10  km  altitude  span  (see  Figure  2.1).  Also  from 
Figures  2.3  (a)  -  2.5  (a)  we  observe  that  the  average  energy  changes  only  by  a  factor  of 
~1.4  (i.e.,  from  ~4.1  to  ~5.7  eV)  with  a  change  in  the  amplitude  of  the  electric  field  of 
a  factor  of  three  (i.e.,  from  5  to  15  V/m).  This  slow  dependence  is  due  to  the  high  value 
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Fig.  2.5.  Results  of  a  single  slab  (at  90 
km  altitude)  simulation  for  a  fixed  ap¬ 
plied  electric  field  of  E  =  5  V/m  applied 
at  time  t  =  0.  Top  left'  Time  evolution 
of  the  normalized  electric  current  and  av¬ 
erage  electron  energy  Top  right  Time 
evolution  of  the  electron  distribution  func¬ 
tion  The  electron  distribution  function  is 
shown  at  0.022  /is,  0.25  /is,  0.75  /is,  2.3 
/is,  14  /is,  and  20  /is.  Bottom  left:  Time 
evolution  of  the  attachment  and  ionization 
rates.  The  current  density  and  excitation 
rates  are  normalized  to  one  electron  per 
cubic  centimeter 


of  the  second  barrier  starting  at  —6  eV  in  the  electon  loss  function  shown  in  Figure  2.2. 
Hence,  the  average  distance  and  time  over  which  the  electrons  travel  in  energy  space 
during  the  heating  by  the  electric  field  starting  from  a  cold  state  can  be  assumed  constant 
for  the  electric  field  amplitudes  discussed  here.  Consequently,  we  conclude  that  over  the 
region  of  strong  attenuation  (80  to  95  km)  of  the  EM  pulses,  the  diffusion  time  for  the 
electron  distribution  (i.e.,  the  time  for  the  distribution  to  reach  equilibrium)  varies  only  by 
a  factor  of  —6  due  to  the  variation  of  the  neutral  density.  We  also  note  that  equilibrium  is 
reached  over  shorter  time  periods  for  the  electric  current,  the  attachment  rate  and  optical 
excitations  with  low  thresholds  than  for  the  ionization  rate  and  optical  excitations  with 
high  thresholds  (see  for  example  Figures  2.3,  2.4,  2.5  (a)  and  (c)  for  comparison).  This 
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fact  allows  us  to  extend  the  validity  of  our  assumption  of  achiving  equilibrium  over  /zs 
time  scales  to  broader  values  of  electric  fields  and  altitudes  as  far  as  the  relationship 
between  the  applied  electric  field  and  the  electric  current  density  is  concerned;  note  that 
the  latter  is  what  is  needed  to  propagate  the  EM  pulses  in  accordance  with  equations 
(2.1).  According  to  single  slab  calculations  at  other  altitudes  not  shown  here,  and  for 
electric  field  amplitudes  of  >5  V/m  that  are  necessary  to  produce  substantial  electron 
density  changes  and  optical  emissions,  equilibrium  is  reached  over  a  several  fis  time 
scale  for  the  current  density,  the  optical  excitation  and  ionization  rates  at  the  altitudes  of 
interest  (80  to  95  km). 

For  the  case  of  small  (<5  V/m)  electric  field  amplitudes,  and  at  altitudes  <85  km, 
the  average  energy  of  the  electrons  in  the  electric  field  is  <2  eV.  For  such  a  case,  the 
most  important  factor  which  determines  the  distribution  is  inelastic  collisions  of  electrons 
with  excitation  of  rotational  levels  and  low  vibrational  levels  (see  for  example  [Gurevich, 
1978;  pp.  70-73]).  The  distribution  function  /  for  such  interactions  is  of  the  Druyvestein 
type  and  decreases  with  velocity  more  rapidly  than  a  Maxwellian  (In  /  ~  IE1).  The 

average  energy  of  such  a  distribution  is  Hence,  the  product  of  the  electric  field 

driving  force  (2.12)  and  the  average  energy  (which  determines  the  time  needed  to  reach 
equilibrium)  varies  as  ~  E 5.  This  analytical  evaluation  shows  that  the  time  needed 
to  reach  equilibrium  decreases  with  electric  field  for  the  lower  electric  field  amplitudes 
and  allows  us  to  state  that  the  fis  time  interval  is  also  sufficient  for  electrons  to  reach 
equilibrium  in  small  electric  fields. 

2.4  MULTIPLE-SLAB  MODEL  OF  THE  INTERACTION  AS  A  FUNCTION 
OF  TIME  AND  ALTITUDE 

The  results  of  our  full  kinetic  solution  of  the  coupling  for  a  single  slab  allow  us  to  formu¬ 
late  a  multi  slab  model  of  the  interaction,  to  describe  the  coupling  and  its  consequences 
as  a  function  of  altitude  and  time,  as  the  EM  pulse  propagates  upward  and  interacts  with 
the  increasingly  dense  lower  ionosphere. 
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Subject  to  the  caveats  discussed  above,  the  major  conclusion  we  derive  from  the  results 
of  the  single  slab  calculations  is  the  fact  that  the  form  of  the  electron  distribution  function 
follows  a  quasi-equilibrium  solution  of  equations  (2.4),  (2.5)  and  is  established  within  a 
time  interval  <20  /x s,  so  that  if  the  applied  electric  field  exhibited  relatively  ‘slow’  field 
changes  over  time  scales  of  tens  of  /xs  then  the  quasi-equilibrium  solution  can  be  used 
to  calculate  the  current  density,  ionization,  and  optical  excitation  rates. 

Available  data  on  the  waveforms  of  lightning  radiated  EM  pulses  (e.g.,  [LJman,  1987; 
p.  113],  also  see  Figure  1.1)  indicate  that  they  exhibit  relatively  slow  time  variations 
over  a  10  /xs  time  scale,  with  pulse  durations  of  typically  50-150  /xs.  Hence,  to  simulate 
the  propagation  of  the  EM  pulses  from  lightning  through  the  lower  ionosphere  we  do 
not  need  to  follow  the  complete  time  evolution  of  the  Boltzmann  equations  (2.4)  and 
(2.5).  The  current  density,  ionization,  attachment,  and  optical  excitation  rates  can  be 
approximated  by  their  values  derived  from  the  equilibrium  solution  for  a  given  electric 
field  existing  at  a  particular  time,  at  a  given  altitude.  For  an  electric  field  increasing 
(decreasing)  with  time  such  an  approximation  leads  to  an  overestimate  (underestimate) 
of  the  number  of  high  energy  electrons  due  to  the  fact  that  the  equilibrium  is  reached 
from  lower  (higher)  energies.  However,  as  we  discuss  below,  the  error  introduced  in  our 
results  due  to  this  approximation  is  typically  less  than  10  %. 

Our  practical  realization  of  the  approach  described  above  was  as  follows.  First,  we 
prepared  a  table  of  equilibrium  values  (calculated  using  the  single  slab  model)  of  current 
density,  electron  density  variation  rates,  and  optical  excitation  rates  as  a  function  of 
altitude  and  electric  field  values.  This  table  covers  the  70-100  km  altitude  range  in  1  km 
increments  and  specifies  the  current  density,  electron  density  variation  rate,  and  optical 
excitation  rates  over  a  50-point  grid  of  values  of  the  electric  field.  The  maximum  value 
of  the  electric  field  over  the  grid  is  70  V/m  below  70  km  altitude  and  is  gradually  reduced 
to  3  V/m  at  100  km.  For  altitudes  above  100  km  the  medium  is  approximated  by  its 
characteristics  at  100  km.  Second,  we  simulated  the  propagation  of  EM  pulses  solving 
Maxwell’s  equations  (2.1)  in  a  self  consistent  manner  choosing  appropriate  values  of 
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the  electric  current,  density  changes,  and  optical  excitation  rates  from  the  table,  using 
linear  interpolation  to  find  values  between  the  tabulated  points  as  shown  in  an  example 
in  Figure  2.6.  Further  details  of  our  numerical  method  are  given  in  Appendix  C. 

Due  to  the  particular  numerical  scheme  employed,  the  present  version  of  the  code 
has  an  upper  limit  for  the  electric  field  amplitude  that  does  not  allow  the  simulation  of 
interactions  with  the  lower  ionosphere  of  EM  pulses  with  amplitude  Einj  >  25.5  V/m 
(observed  ~5%  of  the  time)  for  a  ‘tenuous’  D  region  and  Etnj  >  37.5  V/m  (observed 
~1%  of  the  time)  for  a  ‘dense’  D  region.  Simulation  of  the  interaction  of  such  pulses 
with  the  lower  ionosphere  is  the  subject  of  future  work  (see  Section  6.2). 

To  assess  the  validity  of  our  numerical  technique,  we  monitored  the  electric  field, 
current  density,  and  electron  density  values  at  fixed  altitudes  as  the  EM  pulse  propagated 
through  the  medium  and  then  checked  the  resulting  time  dependences  against  single  slab 
calculations. 

Figure  2.7  (a)  shows  the  variation  of  the  electric  field  at  89  km  calculated  in  the 
manner  described  above  (i.e.,  using  the  grid  values)  for  a  single  cycle  sinusoidal  pulse 
of  100  /is  duration  injected  at  70  km  with  an  initial  peak  intensity  of  Etnj  =  20  V/m  for 
the  ‘tenuous’  nighttime  profile  (a)  of  Figure  2.1.  Figure  2.7  (b)  shows  the  corresponding 
current  values  as  obtained  from  the  same  model  (i.e.,  from  the  corresponding  grid  values), 
and  Figure  2.7  (d)  shows  the  current  density  as  computed  by  a  time  dependent  solution 
of  (2.4),  (2.5)  for  a  single  slab  at  this  altitude  and  for  the  varying  electric  field  value  in 
effect  as  given  in  Figure  2.7  (a).  The  difference  between  the  two  is  barely  visible  on 
the  scale  shown,  confirming  that,  for  the  purposes  of  the  propagation  of  EM  pulses  in 
the  highly  collisional  lower  ionosphere,  our  use  of  the  equilibrium  values  is  a  reasonable 
assumption.  The  changes  in  the  electron  density,  again  calculated  in  the  two  different 
ways,  are  compared  in  Figure  2.7  (c).  The  difference  between  values  calculated  using  the 
equilibrium  solution  and  the  single-slab  time  dependent  solution  with  the  varying  electric 
field  is  ~3.5%.  This  error  is  somewhat  larger  for  smaller  amplitudes  of  the  electric  fields 
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Electric  Field  (V/m) 

Fig.  2.6.  Sketch  of  the  grid  used  to  propagate  EM  pulses  in  the  collisional  ionosphere. 

At  each  point  shown  with  a  dot  the  electric  current  density,  the  rate  of  change  of  electron 
density,  and  the  optical  excitation  rates  were  pre-calculated  using  a  single  slab  model 
similar  to  that  described  in  Section  2.3  at  different  altitudes  and  for  different  applied 
electric  fields.  Linear  interpolation  is  used  between  points  to  find  the  appropriate  values 
of  quantities  in  intermediate  points  (e.g.,  data  from  points  1  to  4  are  used  to  perform 
calculations  at  the  point  designated  with  +).  In  the  actual  grid  used  for  our  model,  the 
grid  points  are  separated  by  1  km  in  altitude  (i.e.,  30  points  between  70  and  100  km)  and 
the  electric  field  was  represented  over  a  50-point  grid,  with  the  maximum  value  varying 
from  70  V/m  at  70  km  to  3  V/m  at  100  km. 


since  the  system  reaches  equilibrium  over  a  longer  time  intervals  for  weaker  driving  fields. 
However,  since  for  injected  electric  field  amplitudes  E{n)  <  10  V/m  the  resultant  electron 
density  changes  are  <1  %  of  the  ambient  we  can  safely  use  the  approach  described  above 
for  the  range  of  parameters  of  interest.  We  note  that  during  the  calculations  the  ratio  of 
the  average  directional  electron  velocity  to  the  random  electron  velocity  was  0.01  to  0.1 
which  confirms  the  validity  of  expansion  (2.3). 

Overall,  the  error  introduced  in  our  estimates  of  electron  density  changes  because 
of  our  usage  of  the  quasi-equilibrium  solution  and  the  method  described  above  is  <10 
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Time,  psec  Time,  psec 

Fig.  2.7.  Dynamic  test  of  the  numerical  technique.  Variations  with  time  of  the  electric 
field  (a)  and  the  electric  current  density  (b)  at  89  km  altitude  for  the  ‘tenuous’  nighttime 
conditions  (profile  (a)  of  Figure  2.1)  following  the  injection  of  a  Ein]  =  20  V/m  EM 
pulse  starting  at  t  =  0  at  70  km.  Single  slab  calculations  of  the  electric  current  for  the 
electric  field  varying  as  shown  in  (a)  at  89  km  are  shown  in  plot  (d).  Plot  (c)  shows  a 
comparison  of  the  electron  density  changes  as  calculated  from  the  grid  model  (solid  line) 
and  from  the  single  slab  calculations  (lighter  line)  at  89  km  for  the  electric  field  varying 
as  shown  in  (a). 


%.  Due  to  a  monotonic  decrease  in  the  electron  distribution  function  with  energy  (see 
e.g..  Figure  2.3  (b))  and  the  similarity  of  the  single  peaked  dependences  of  the  electron 
impact  excitation  cross  sections  of  optical  electronic  states  and  dissociative  attachment 
and  ionization  processes,  the  error  in  optical  excitation  rates  would  be  similar  (<10  %). 
Such  accuracy  is  quite  acceptable  for  the  problem  at  hand,  especially  in  view  of  the  fact 
that  the  absolute  magnitudes  of  the  input  radiation  fields  from  lightning  exhibit  a  much 
wider  dynamic  range,  and  that  the  ambient  nighttime  D  region  electron  density  is  also 
highly  variable  and  not  well  known  [Friedrich  and  Torkar,  1992;  Mitra,  1981]. 
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2.5  REPRESENTATION  OF  INJECTED  EM  PULSES  FROM  LIGHTNING 

To  put  the  input  electric  field  amplitudes  E,nj  in  the  right  context,  it  suffices  to  remember 
that  the  mean  amplitude  for  CG  flashes  as  observed  at  100  km  distance  on  the  ground 
is  Ejoo  ~8.8  to  11.2  V/m,  with  amplitudes  >20  V/m  occuring  ~3  to  8  %  of  the  time 
[Krider  and  Guo,  1983]. 

Since  there  is  no  significant  attenuation  below  80  km  for  parameters  of  interest  (Chap¬ 
ter  3)  and  since  we  use  a  plane  wave  model  in  our  calculations,  it  is  useful  to  refer 
the  injected  EM  pulse  amplitude  ( Einj )  to  this  altitude.  The  lack  of  spherical  spreading 
losses  due  to  the  plane  wave  assumption  introduces  ~10  %  error  in  our  calculations, 
which  is  comparable  with  other  uncertainties  due  to  other  assumptions  introduced  in  the 
model,  especially  in  view  of  the  fact  that  the  absolute  magnitudes  of  the  input  radiation 
fields  from  lightning  exhibit  a  much  wider  dynamic  range  and  that  practically  all  the 
data  on  the  lightning  electric  field  amplitudes  come  from  the  ground  based  measurements 
[Uman,  1987;  p.  110].  We  note  that  there  is  reason  to  believe  that  lightning  radiation 
fields  above  the  ground  may  be  substantially  different  [Krider,  1992]. 

Assuming  a  spherically  spreading  wave,  and  leaving  aside  directional  aspects  of  the 
radiation  from  CG  versus  IC  discharges  [e.g.,  see  Rodriguez  et  al.,  1992]  (i.e.,  assuming 
an  isotropic  radiator),  an  electric  field  amplitude  of  Eioo  would  be  related  to  an  injected 
field  amplitude  EtrtJ  at  80  km  as  Etnj  x  ~|g.  In  accordance  with  this,  the  plane  EM 
wave  can  be  injected  at  any  altitude  below  80  km  with  the  same  amplitude  E{n}  without 
significant  error  (note  that  there  is  no  spherical  spreading  factor  in  our  model);  for  results 
presented  in  this  thesis,  we  adopt  an  injection  altitude  of  70  km.  We  emphasize  that 
moving  the  injection  altitude  to  any  other  altitude  below  70  km  does  not  change  any  of 
the  results  of  our  simulations  due  to  the  plane  wave  model  and  lack  of  attenuation  at 
those  altitudes. 

The  optical  emission  intensities  and  the  altitude  distribution  and  absolute  magnitude  of 
ionization  changes  are  highly  sensitive  to  the  time  duration  of  the  lightning  EM  pulse,  its 
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shape  and  intensity.  Each  of  these  lightning  characteristics  exhibits  great  variability  for 
both  CG  and  IC  discharges;  see  examples  in  Figure  1.1.  For  example,  the  time  duration 
of  the  first  return  stroke  varies  by  an  order  of  magnitude  with  ~100  //s  being  the  average 
[Uman,  1987;  p.  110],  and  the  amplitude  of  the  radiated  EM  field  can  reach  E\qo  >25 
V/m  with  the  mean  value  -Ejoo  ~10  V/m  [ Krider  and  Guo ,  1983].  The  shape  of  the 
pulse  is  also  highly  variable,  including  both  single  and  double  polarities,  which  produce 
different  extremes  of  the  maximum  electric  field  amplitude  in  the  lower  ionosphere  when 
the  upgoing  tail  of  the  pulse  interferes  with  the  ionospherically  reflected  front  part  of 
the  pulse.  In  this  thesis  our  objective  is  to  accurately  determine  the  intensities  of  optical 
emissions  and  electron  density  changes  for  EM  pulses  with  a  time  duration  and  shape 
for  typical  lightning  waveform;  investigation  of  the  coupling  of  the  full  range  of  possible 
waveforms  is  the  subject  of  future  work  (see  Section  6.2).  For  this  purpose,  we  represent 
the  injected  EM  pulses  by  a  single  period  oscillation  of  a  10  kHz  sinusoid.  The  largest 
electric  field  amplitudes  ( Ein]  =  25.5  V/m  (observed  ~5%  of  the  time)  for  profile  (a) 
and  Etnj  =  37.5  V/m  (observed  ~1%  of  the  time)  for  profile  (b))  that  are  used  in  our 
calculations  are  limited  by  the  particular  numerical  model  employed. 


3 

Electron  Density  Changes 


To  study  the  effect  on  the  ionosphere  of  radiation  from  a  single  lightning  discharge  we 
inject  a  100  /is  long  EM  pulse  at  the  lower  boundary  (70  km  altitude)  represented  by 
a  single  period  oscillation  of  a  10  kHz  sinusoidal  wave.  We  simulate  the  propagation 
of  the  EM  pulse  solving  Maxwell’s  equations  (2.1)  and  its  evolution  as  it  interacts  with 
the  ionosphere  as  determined  by  equations  (2.2)-(2.6)  for  a  time  duration  of  ~260  //s, 
which  allows  the  injected  pulse  to  be  reflected  from  the  lower  ionospheric  boundary  (~90 
km)  and  to  propagate  back  through  the  interaction  region  on  its  way  to  the  ground.  We 
maintain  zero  field  boundary  conditions*  at  the  upper  boundary  (120  km)  consistent  with 
strong  attenuation  and  reflection  in  the  lower  layers.  The  resulting  density  perturbations 
are  shown  at  the  end  of  the  interaction  for  a  variety  of  field  amplitudes  and  ambient 
electron  density  profiles.  We  also  illustrate  a  time  sequence  of  events  with  snap-shots 
of  altitude  distributions  of  electric  field,  electron  density  changes,  and  average  energy 
of  electrons  for  the  interaction  of  a  Exnj  =  20  V/m  pulse  with  the  ‘tenuous’  nighttime 
ionosphere. 

Often  the  time  period  between  consecutive  discharges  from  thunderstorm  centers  is  less 

*  We  also  noie  that  due  to  the  quasi-steady  ( d/dt  ~  0)  approach  to  the  solution  of  the  Boltzmann  equations 
(2.4)-{2.6)  our  model  does  not  account  for  whistler-mode  waves  (for  any  orientation  of  the  electric  field) 
which  penetrate  to  the  upper  ionosphere  and  propagate  to  the  weakly  collisional  magnetized  plasma 
of  the  magnetosphere  [e.g.,  Helliwell,  1965].  We  note  however,  that  leaving  the  whistler-mode  out  of 
our  model  has  no  significant  impact  on  our  modeling  of  the  interaction  of  EM  pulses  with  the  highly 
collisional  plasma  of  the  lower  ionosphere. 
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than  a  second  [Uman,  1987,  p.  19],  which  is  much  less  than  the  relaxation  time  (10  to  100 
s)  of  the  density  perturbations  at  D-region  altitudes  [Glukhov  et  al.,  1992].  In  such  cases, 
the  ionization  produced  by  consecutive  strokes  would  accumulate  as  previously  suggested 
[Inan  et  al.,  1991].  In  modeling  this  accumulation  we  neglect  relaxation  during  the  time 
between  strokes  and  supply  incident  EM  energy  for  ~2.2  ms,  a  time  that  corresponds  to 
~20  consecutive  lightning  strokes.  The  resulting  density  perturbations  are  also  presented 
for  different  electron  densities  and  field  strengths  ( E  ||  Bo).  First  we  present  the  results 
for  the  electric  field  of  the  pulse  oriented  parallel  to  the  ambient  magnetic  field.  At  the 
end  of  the  chapter  we  compare  these  results  with  results  for  the  perpendicular  orientation 
(E  ±  B0). 

3.1  ELECTRON  DENSITY  CHANGES  FOR  A  ‘TENUOUS’  AMBIENT 
PROFILE 

We  start  with  the  ambient  electron  density  profile  (a)  in  Figure  2.1  corresponding  to  a 
relatively  ‘tenuous’  nighttime  ionosphere.  Figure  3.1  shows  snapshots  covering  the  80  to 
180  /xs  time  interval  in  20  /xs  increments  following  the  injection  of  the  EM  pulse  showing 
the  altitude  distribution  of  the  electric  field,  average  electron  energy,  and  normalized 
electron  density  perturbations  (SNe/Ne o). 

The  injected  wave  amplitude  at  70  km  was  taken  to  be  Einj  =  20  V/m  with  the  electric 
field  polarized  parallel  to  the  ambient  magnetic  field  (assumed  horizontal).  At  80  /xs  the 
wave  front  reaches  94  km  altitude.  As  the  wave  propagates  it  becomes  progressively 
attenuated  and  is  partially  reflected.  We  see  from  Figure  3.1  at  80  /xs  that  at  85  km  the 
amplitude  has  decreased  to  ~  17.5  V/m.  However,  due  to  the  partial  reflection  from  the 
lower  ionospheric  boundary  the  amplitude  of  the  pulse  in  the  medium  can  exceed  the 
initial  amplitude  ( E{nj  =  20  V/m  in  this  case)  as  a  result  of  constructive  interference. 
An  example  of  this  effect  is  seen  at  120  /xs  as  the  amplitude  reaches  ~27  V/m  at  ~83 
km.  Examination  of  Figure  3.1  shows  that  the  lightning  EM  fields  which  penetrate  above 
97  km  are  strongly  attenuated.  The  average  energy  of  the  electrons  at  ~90  km  altitudes 
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>.  3.1.  Dynamics  of  the  interaction  of  an  EM  pulse  with  the  ionosphere.  Snapshots 
the  electric  field,  normalized  electron  density  change,  and  average  electron  energy 
Bering  the  80  to  180  /is  time  interval  in  20  ps  increments  following  the  injection  of 
5i„j  =  20  V/m  EM  pulse  at  70  km  under  the  ‘tenuous’  nighttime  conditions  (Fig- 
:  2.1  (a)). 
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at  times  exceeds  5  eV,  providing  a  substantial  number  of  electrons  with  energy  higher 
than  the  thresholds  for  dissociative  attachment  (~6  eV)  and  ionization  (~16  eV)  and 
thus  leading  to  changes  in  electron  density.  At  lower  altitudes  (79  to  86  km)  and  at  the 
end  of  the  region  penetrated  by  the  fields  (92  to  95  km)  dissociative  attachment  prevails, 
causing  decreases  in  electron  density.  At  intermediate  altitudes  (86  to  92  km)  the  density 
is  increased  as  ionization  dominates.  The  major  increase  in  density  occurs  at  the  time 
of  constructive  interference  of  the  downgoing  reflected  first  half  of  the  pulse  with  the 
upgoing  second  half.  The  value  of  the  density  change  reaches  ~-3  %  at  83  km  and  >+7 
%  at  89  km  at  the  end  of  the  interaction. 

The  resulting  density  changes  for  6,  12,  and  18  successive  EM  pulses  with  Etn]  =  20 
V/m  are  shown  in  Figure  3.2.  The  depletion  evolves  to  ~48  %  of  the  ambient  at  83  km 
and  the  accumulated  density  increase  is  >150  %  at  89  km. 

A  single  100  /z s  pulse  with  E{nj  =  25  V/m  causes  up  to  a  27  %  increase  at  ~89  km 
and  a  decrease  of  ~4  %  at  ~83  km  as  shown  in  Figure  3.3  (a).  The  resulting  density 
changes  for  8,  14,  and  20  successive  lightning  strokes  with  Etn]  =  25  V/m  are  shown 
in  Figure  3.3  (b).  The  depletion  evolves  to  ~55  %  of  the  ambient  at  ~83  km  and  the 
density  increase  is  >350  %  at  89  km. 

A  summary  plot  of  the  resulting  normalized  density  changes  for  single  pulses  of  differ¬ 
ent  strength  is  shown  in  Figure  3.4.  Note  for  example  that  a  pulse  with  100  /xs  duration 
but  with  a  smaller  initial  amplitude  E{nj= 10  V/m  causes  a  substantially  smaller  (<1  %) 
negative  density  change  with  its  maximum  located  at  ~87  km,  whereas  a  ~20  %  density 
decrease  occurs  for  20  successive  EM  pulses  of  this  strength.  The  drastic  differences  be¬ 
tween  the  results  for  initial  field  intensities  of  10,  20,  and  25  V/m  underscore  the  highly 
nonlinear  nature  of  the  interaction. 

It  is  also  interesting  to  note  from  the  left  hand  panel  of  Figure  3.4that  the  altitude 
of  maximum  depletion  moves  lower  with  increasing  input  electric  field  intensity.  This 
tendency  is  opposite  to  that  found  in  earlier  work,  namely  that  the  altitude  of  maximum 
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Fig.  3.2.  Electron  density  changes.  The  resulting  density  changes  produced  by  a  se¬ 
quence  of  8,  14,  and  18  EM  pulses  with  £i„;  =  20  V/m  under  ‘tenuous’  nighttime 
conditions  (Figure  2.1  (a)). 


heating  increases  with  injected  EM  pulse  intensity  [Rodriguez  et  al.,  1992],  due  to  the 
increase  in  reflection  height  as  a  result  of  enhanced  collision  frequency  (in  the  context 
of  a  WKB  approximation  the  reflection  of  VLF  waves  occurs  at  the  altitude  at  which 
~  uv  [Ratcliffe,  1959;  p.  Ill]  where  up  is  the  plasma  frequency,  u  is  the  wave 
frequency,  and  v  is  the  effective  electron  collision  frequency). 
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Fig.  3.3.  Electron  density  changes,  (a):  The  resulting  density  changes  produced  by 
a  single  EM  pulse  with  Einj  =  25  V/m  under  ‘tenuous’  nighttime  conditions  (Fig¬ 
ure  2.1  (a)),  (b):  The  resulting  density  changes  produced  by  a  sequence  of  8,  14,  and  20 
EM  pulses  with  £,n;  =  25  V/m  under  the  same  conditions. 


3.2  ELECTRON  DENSITY  CHANGES  FOR  A  ‘DENSE’  AMBIENT  PROFILE 

As  an  intermediate  ambient  density  profile,  we  consider  profile  (b)  in  Figure  2.1,  which 
may  represent  nighttime  conditions  during  enhanced  geomagnetic  activity  at  midlatitude 
or  in  the  auroral  regions. 

For  a  denser  D  region  the  reflection  altitude  for  the  EM  pulse  is  lower,  resulting  in 
higher  local  neutral  density  and  less  electron  heating  compared  to  a  tenuous  D  region. 
Thus,  larger  electric  field  amplitudes  are  needed  to  produce  the  same  fractional  density 
charges  or  to  excite  a  given  intensity  of  optical  emissions. 

For  an  EM  pulse  with  E{nj  =  35  V/m  the  density  changes  are  shown  in  Figure  3.5. 
Negative  changes  (~-4  %  at  81  km)  are  of  larger  magnitude  than  positive  ones  (~+2.3 
%  at  87  km).  For  pulses  with  Exn]  <30  V/m,  the  density  perturbations  become  single 
peaked  and  purely  negative  at  ~83  km.  The  absolute  value  of  8Ne/Ne  drops  below  1  % 
for  EM  pulses  with  initial  amplitudes  smaller  than  20  V/m. 


Altitude  (km) 
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Fig.  3.4.  Electron  density  changes.  Summary  of  the  resulting  density  changes  produced 
by  single  EM  pulses  of  Etnj  =  10, 12.5, 15, 17.5,20,22.5,25,25.5  V/m  under  ‘tenuous’ 
nighttime  conditions  (Figure  2.1  (a)).  Note  the  different  scales  for  the  abscissa  for  the 
left  and  right  hand  panels. 

Successive  EM  pulses  with  initial  amplitude  Exnj  =  35  V/m  can  produce  quite  promi¬ 
nent  density  changes  as  shown  in  Figure  3.6  for  8,  14,  and  20  discharges.  The  positive 
and  negative  changes  have  about  the  same  amplitude  (~60  %)  and  peak  at  85  km  and 
81  km  respectively.  For  initial  electric  field  amplitudes  Etnj  =30  V/m,  the  density  per¬ 
turbation  is  negative  with  a  single  peak  at  about  83  km  and  an  amplitude  of  ~57  %  for 
Einj  =  30  V/m,  and  ~20  %  for  Exn]  -  20  V/m. 

3.3  DAYTIME  CONDITIONS 

The  daytime  profile  shown  in  Figure  2.1  as  profile  (c)  was  used  to  make  calculations 
similar  to  those  reported  above.  In  general,  we  find  that  even  for  the  strongest  EM 
pulses  ( E{nj  =  50  V/m)  the  density  perturbation  is  negative  with  very  small  absolute 
value  (\6Ne/Ne\  <  3  x  10-?),  whereas  for  20  consecutive  discharges  of  such  intensity 
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Fig.  3.5.  Electron  density  changes.  The  resulting  density  changes  produced  by  a  single 
Einj  =  35  V/m  EM  pulse  for  a  ‘dense’  electron  density  profile  (Figure  2.1  (b)). 

( 1 6 Ne /Ne  |  <  5  x  10-2).  The  minimum  is  located  at  ~76  km.  Such  perturbations 
are  unlikely  to  be  detectable  and  most  probably  do  not  play  a  significant  role  in  the 
overall  dynamics  of  the  lower  ionosphere.  We  therefore  conclude  that  the  phenomenon 
of  lightning-induced  ionization  occurs  predominantly  under  nighttime  conditions. 

The  physical  reason  for  the  lack  of  significant  ionization  by  lightning  during  daytime 
is  the  fact  that  the  reflection  altitude  of  the  lightning  EM  pulses  is  at  ~77  km,  where 
the  mean  free  path  of  the  electrons  is  much  shorter  (~4  times)  than  at  85  km,  and 
consequently,  more  of  the  electron  energy  goes  into  excitation  of  the  molecular  energy 
levels  lower  than  those  of  ionization  and  dissociative  attachment. 
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Fig.  3.6.  Electron  density  changes.  The  resulting  density  changes  produced  by  a  se¬ 
quence  of  8,  14,  and  20  Etnj  =  35  V/m  EM  pulses  under  the  ‘dense’  nighttime  profile 
(Figure  2.1  (a)). 

3.4  DEPENDENCE  ON  THE  ORIENTATION  OF  THE  ELECTRIC  FIELD 

The  results  reported  so  far  were  for  an  EM  pulse  with  electric  field  assumed  to  be 
parallel  to  the  ambient  magnetic  field  ( E  _L  Dq),  which  in  turn  is  taken  to  be  horizontal. 
After  performing  similar  calculations  for  the  case  of  the  perpendicular  orientation  of 
the  electric  field  of  the  pulse  with  respect  to  the  ambient  magnetic  field,  we  found  two 
minor  differences  from  the  case  of  the  parallel  field  orientation.  Firstly ,  for  the  same 
electric  field  amplitudes  and  pulse  duration,  the  maximum  amplitudes  of  the  positive 
density  perturbations  for  the  perpendicular  orientation  were  15  to  30  %  smaller  than  for 
the  case  of  parallel  orientation.  Secondly,  the  altitudes  of  the  maximum  positive  density 
perturbations  for  the  perpendicular  orientation  were  2  to  3  km  higher  than  for  the  case  of 
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parallel  orientation.  For  example,  a  E{nj  =  20  V/m  100  /x s  pulse  with  the  perpendicular 
orientation  for  the  nighttime  conditions  produced  ~6  %  density  increase  located  at  ~91 
km  in  comparison  to  >7  %  increase  at  89  km  for  the  parallel  orientation  (see  Figure  3.1). 

This  comparatively  weak  dependence  on  the  angle  between  the  electric  field  of  the 
wave  and  the  ambient  magnetic  field  is  due  to  the  high  collision  frequency  of  the  heated 
electrons  with  the  ambient  neutrals  which  is  much  larger  than  the  electron  gyrofrequency 
in  the  ambient  magnetic  field  in  the  altitude  range  of  interest.  For  example,  at  90  km  the 
collision  frequency  of  5  eV  electrons  is  ~  108  s-1  which  is  about  two  orders  of  magnitude 
larger  than  the  electron  gyrofrequency  at  that  altitude  which  makes  the  influence  of  the 
magnetic  field  completely  negligible.  However,  for  weaker  heating  (for  example  in  the 
case  of  heating  by  VLF  transmitters)  D0  is  important  and  leads  to  substantial  asymmetry 
of  the  transverse  shape  of  the  heated  region  [Inan  et  al,  1992;  Rodriguez,  1993], 

Overall,  the  results  obtained  in  this  study  appear  to  be  applicable  with  ~10  %  accuracy 
to  any  intermediate  (between  E  ||  Bq  and  E  _L  Bq)  field  orientation. 


4 

Optical  Emissions 


Airglow  and  auroral  emissions  have  been  used  extensively  to  extract  information  about 
important  ionospheric  and  magnetospheric  processes.  Excitation  of  a  variety  of  optical 
emissions  would  be  a  natural  consequence  of  lightning-induced  heating  of  ambient  elec¬ 
trons  as  suggested  by  lnan  et  al.  [1991]  and  as  discussed  in  Chapter  3.  The  experimental 
data  that  motivated  the  work  of  lnan  et  al.  [1991  ]  involved  detection  (via  subionospheric 
VLF  remote  sensing)  of  lightning-associated  ionization  enhancements  (as  discussed  in 
Chapter  3)  that  are  produced  at  threshold  energies  of  >  15.6  eV,  which  is  well  above  the 
excitation  thresholds  of  most  optical  emissions  from  neutral  constitutients  in  the  D  region 
[e.g.,  Rees,  1989;  Appendix  3].  The  possibility  of  the  excitation  of  optical  emissions  as  a 
result  of  lightning-induced  heating  of  the  lower  ionosphere  was  first  recognized  by  lnan 
et  al.  [1991],  The  first  experimental  evidence  for  such  optical  emissions  via  D  region 
heating  by  lightning  may  be  the  recent  observations  from  the  Space  Shuttle  [Boeck  et 
al.,  1992]  of  a  transient  airglow  enhancement  in  the  D  region  of  the  ionosphere  that 
appeared  to  be  associated  with  a  tropospheric  lightning  discharge.  These  observations 
were  discussed  earlier  in  Section  1 .2. 

In  this  Chapter  we  use  the  same  model  for  the  electrodynamic  interaction  of  a  lightning 
EM  pulse  with  the  lower  ionosphere  as  that  used  in  Chapter  3  to  calculate  the  radiation 
time  constants  and  intensities  of  the  emission  lines  of  O,  O2,  N2,  Oj,  and  N}  most  easily 
detectable  in  the  natural  aurora  [Vallance  Jones,  1974], 
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To  study  optical  emissions  from  the  D  region  excited  by  EM  radiation  from  a  single 
lightning  stroke  we  inject,  as  before  in  Chapter  3,  a  100  /zs  long  EM  pulse  at  the  lower 
boundary  (70  km  altitude)  represented  by  a  single  period  oscillation  of  a  10  kHz  sinusoidal 
wave.  In  accordance  with  the  model  described  in  Chapter  2,  we  simulate  the  propagation 
of  the  pulse  solving  Maxwell’s  equations  (2.1)  and  its  evolution  as  it  interacts  with 
the  ionosphere  as  determined  by  equations  (2.2)-(2.6)  for  a  time  duration  of  ~260  /zs, 
which  allows  the  injected  pulse  to  be  reflected  from  the  lower  ionospheric  boundary 
(~90  km)  and  to  propagate  back  through  the  interaction  region  on  its  way  to  the  ground. 
We  maintain  zero  field  conditions  at  the  upper  boundary  (120  km)  consistent  with  strong 
attenuation  and  reflection  in  the  lower  layers.  In  this  Chapter  we  present  time  snapshots  of 
the  altitude  distributions  of  electric  fields,  average  electron  energy  and  emission  intensity, 
altitude  distributions  of  the  maximum  intensities  of  optical  emissions,  and  time  and  height 
integrated  emission  intensities. 

The  governing  equation  in  our  model  for  the  impact  excitation  of  neutral  species  is 
[e.g.,  Sipler  and  Biondi,  1972] 

^r  =  -^+T,n>Ai+N-R-  '4I> 

*  j 

where  n,  denotes  the  density  of  excited  particles  at  z-th  level,  N,  is  the  ambient  density 
of  corresponding  neutrals  in  the  ground  state,  Ti  =  (£^  Ai+kuNN  +k2iN02+k3,N0)~l  is 
the  total  lifetime  of  the  z-th  state,  A\  is  the  radiation  transition  rate*,  kij  is  the  quenching 
rate!,  and  R,  is  the  excited  state  source  term4  The  first  term  on  the  right  hand  side  of 
equation  (4.1)  describes  the  decrease  of  the  density  of  the  excited  particles  as  a  result 
of  emissions  and  quenching  by  neutrals,  the  second  term  represents  the  increase  in  the 
population  of  excited  particles  as  a  result  of  cascading  from  higher  states  (e.g.,  population 

*  The  transition  rate  is  the  inverse  of  the  average  time  it  takes  for  the  spontaneous  radiation  of  a  photon. 

t  Quenching  refers  to  the  loss  of  energy  by  the  excited  particle  due  to  collision  processes  with  molecules 
or  atoms. 

$  The  source  term  describes  excitation  of  the  state  by  electron  impact  from  the  ground  state. 
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of  the  B^ng  state  of  N2  as  a  result  of  emissions  of  the  2nd  positive  band  originating  from 
the  C37ru  level  in  Figure  4.1),  and  the  third  term  describes  the  excitation  of  a  particular 
state  by  electron  impact  from  the  ground  state  (e.g.,  X'l^  state  of  N2  in  Figure  4.1). 


Fig.  4.1.  Energy  levels  and  electronic  band  systems  for  N2  and  NJ.  The  1st  and  2nd 
positive  bands  of  N2  and  1st  negative  and  Meinel  bands  of  Nj  are  considered  in  our 
calculations.  For  a  variety  of  different  reasons,  other  bands  of  N2  and  Nj  (not  shown 
here)  are  considered  to  be  weak  based  on  observed  auroral  spectra  [Vallance  Jones,  1974; 
p.  90], 


The  intensity  of  t-th  line  is  given  by  [Chamberlain,  1978;  p.  213]: 

/,  =  1(T6  J  Atnxdz  (4.2) 

where  the  integration  is  performed  along  the  line  of  sight  and  the  units  of  the  resultant 
intensity  /,  are  in  Rayleighs  (R),  (i.e.,  photons/s  emitted  in  a  column  with  1  cm2  cross 
sectional  area  along  the  line  of  sight  multiplied  by  10-6). 
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Differences  between  the  time  duration  of  EM  pulses  (~100  /xs)  and  the  time  intervals 
of  spontaneous  radiation  and  quenching  allow  us  to  consider  the  process  of  excitation  in 
a  relatively  simple  manner.  The  fact  that  the  100  /rs  time  interval  is  much  shorter  than 
the  emission  and  quenching  times  for  the  ‘forbidden’*  [Vallance  Jones,  1974;  p.  95] 
transitions  (e.g.,  6300,  5577  A  of  O  with  the  quenching  times  being  1.5  ms  and  0.35 
s  at  ~95  km,  emission  times  being  107  s  and  1  s  (Appendix  D)  respectively),  allows 
us  to  neglect  emission  and  quenching  losses  during  the  time  interval  of  the  excitation 
process  and  to  integrate  the  number  of  excited  molecules  over  the  time  of  excitation.  The 
resulting  number  of  excited  atoms  is  used  to  calculate  the  maximum  intensities  of  the 
lines.  For  all  molecular  bands  (which  constitute  ‘allowed’t  [Vallance  Jones,  1974;  p.  95] 
transitions),  the  fact  that  the  100  /is  time  interval  is  much  longer  than  the  time  intervals 
(the  longest  is  14  /is  for  Meinel  band  of  N2,  see  Appendix  D)  of  spontaneous  radiation  of 
these  emissions  and  much  shorter  than  the  quenching  time  intervals  ~1  ms  (Appendix  D) 
at  altitudes  of  interest  allows  us  to  assume  that  ~100  %  of  the  excited  molecules  emit 
photons  at  about  the  excitation  rates  in  effect  during  the  time  of  the  heating. 

The  three  major  constituents  for  which  we  estimate  the  optical  emission  characteristics 
are  N2,  O2,  and  O.  We  also  consider  O2  and  N2  ions  produced  as  a  result  of  the  electron 
impact  ionization  of  the  neutral  species.  Energy  level  diagrams  of  O,  O2,  N2,  O2,  and 
are  well  known  and  are  widely  available  in  literature  (e.g.,  see  Vallance  Jones  [1974], 
pp.  85-91).  The  N2  energy  diagram  is  given  as  an  example  in  Figure  4.1. 

Each  of  the  neutral  and  ionized  species  under  consideration  has  emission  lines  ranging 
from  far  infrared  (IR)  to  extreme  ultraviolet  (UV)  which  are  commonly  observed  in 
auroral  and  airglow  emission  spectra.  Emissions  to  be  studied  here  were  selected  upon 


*  The  term  ‘forbidden’  is  used  to  indicate  the  fact  that  these  states  have  a  zero  electric  and/or  magnetic 
dipole  moment  as  a  consequence  of  the  symmetry  properties  of  the  quantum  wave  function. 

t  The  term  ‘allowed’  is  used  to  indicate  the  fact  that  these  states  have  a  electric  or  magnetic  dipole  moment 
not  equal  to  zero  as  a  consequence  of  the  symmetry  properties  of  the  quantum  wave  function. 
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consideration  of  the  following  factors:  i)  amplitudes  of  electron  impact  excitation  cross- 
sections  (electron  impact  transfers  kinetic  energy  of  the  electrons  to  optical  excitations 
of  ambient  molecules;  the  individual  cross  sections  are  summarized  in  Appendix  B),  ii) 
emission  rates  (which  determine  the  average  time  intervals  of  spontaneuous  emissions 
summarized  in  Appendix  D),  iii)  branching  ratio  of  different  lines  originating  from  the 
same  level  (e.g.,  the  ratio  of  the  number  of  5577  A  ('S  — AD  transition)  and  2972  A 
('5  — transition)  photons  originating  from  the  same  'S  level  of  O;  see  Appendix  D), 

iv)  cascade  effects  (e.g.,  population  of  the  B*ng  state  of  N2  as  a  result  of  emissions  of 
the  2nd  positive  band  originating  from  the  C3 nu  level  as  can  be  seen  from  Figure  4.1), 

v)  radiation  trapping  effects  (i.e.,  resonant  absorption  of  emitted  photons  by  the  same 
type  of  atoms  or  molecules,  e.g.,  1304  and  1356  A  lines  of  O  [Donahue  and  Strickland, 
1970]),  vi)  quenching  rates  by  neutral  species  (upon  collision  of  an  electronically  excited 
molecule  with  another  molecule  (or  atom)  in  a  ground  state  the  former  can  lose  its 
energy  to  the  excitation  of  a  vibrational  state  or  to  the  translational  energy  of  the  latter 
[Vallance  Jones,  1974;  p.  114];  the  quenching  rates  are  summarized  in  Appendix  D), 
vii)  transparency  of  the  atmosphere  at  the  wavelength  of  the  emission  (after  elimination 
of  lines  resonantly  absorbed  in  the  atmosphere,  the  major  factor  determining  the  optical 
depth  of  the  path  is  the  Rayleigh  molecular  scattering  discussed  in  Section  5.2.2),  viii) 
comparison  of  the  expected  intensities  to  the  nighttime  airglow  background  (this  aspect 
is  discussed  in  Section  5.2.1). 

The  optical  lines  and  bands  considered  here  are  regularly  observed  in  the  airglow  and 
auroral  spectra  [Vallance  Jones,  1974;  Chamberlain,  1978],  namely:  the  red  6300  (6364) 
\(]D  — >  3P  transitions),  and  the  green  5577  A  ('  5  — >  1 D  transitions)  lines  of  O;  the  blue 
to  UV  2nd  positive  band  (C37t„  — ♦  B^Xg  transitions)  and  the  red  to  IR  1st  positive  band 
(B3 -Kg  — *  A3I+  transitions)  of  N2;  the  blue  1st  negative  band  (B2£+  — >  X2lAg  transition) 
and  IR  Meinel  band  (A2ttu  — >  X21^  transition)  of  N2  ions;  and  the  red  1st  negative 
band  (64Z“  — *  a4  7ru)  of  Oj  ions.  We  exclude  the  Herzberg  I  band  from  our  present 
analysis  due  to  lack  of  reliable  electron  impact  excitation  cross  sections  for  this  band  [T. 
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Slanger,  private  communication].  The  O2  atmospheric  band  is  also  excluded  due  to  the 
weak  emission  intensities  in  comparison  with  the  nighttime  backgrounds  [Taranenko  et 
al.,  1992].  Plots  of  the  cross-sections  for  excitation  of  the  chosen  emissions  by  electron 
impact  as  a  function  of  electron  energy  are  shown  in  Appendix  B  whereas  quenching  and 
emission  rates  are  summarized  in  Appendix  D. 

4.1  RESULTS  FOR  A  ‘TENUOUS’  NIGHTTIME  D  REGION 

We  start  with  the  ambient  electron  density  profile  (a)  in  Figure  2.1  corresponding  to  a 
‘tenuous’  nighttime  D  region  and  show  the  dynamics  of  the  interaction  of  the  EM  pulse 
with  the  lower  ionosphere. 

Figure  4.2  shows  snapshots  covering  the  80  to  180  [is  time  interval  in  20  fis  increments 
following  the  injection  of  the  EM  pulse.  Each  frame  shows  the  altitude  distribution  of 
the  electric  field,  average  electron  energy,  and  number  of  emissions  per  cc  per  second  of 
the  5577  A  line  of  O  and  the  1st  positive  band  of  N2.  The  left  and  right  hand  panels 
of  Figure  4.2  are  identical  to  those  of  Figure  3.1;  the  center  panels  in  Figure  4.2  show 
optical  emission  intensities  instead  of  density  changes  as  in  Figure  3.1. 

The  injected  wave  amplitude  at  70  km  was  taken  to  be  20  V/m  with  the  electric 
field  polarized  parallel  to  the  ambient  magnetic  field  (assumed  horizontal).  At  80  yas 
the  wave  front  reaches  94  km  altitude.  As  the  wave  propagates  upward  it  becomes 
progressively  attenuated  and  is  partially  reflected.  The  average  energy  of  the  electrons  at 
times  exceeds  5  eV,  providing  a  substantial  number  of  electrons  with  energy  higher  than 
the  optical  thresholds  (ranging  from  2  to  19  eV)  and  thus  leading  to  stimulation  of  optical 
emissions.  It  is  seen  that  the  intensity  of  the  ‘forbidden’  line  (5577  A)  grows  throughout 
the  process  of  interaction  whereas  the  intensity  of  the  1st  positive  ‘allowed’  band  is  highly 
variable  and  reaches  its  largest  amplitude  at  the  time  of  constructive  interference  of  the 
downgoing  reflected  first  half  of  the  pulse  with  the  upgoing  second  half.  The  difference 
in  their  behaviour  is  due  to  the  difference  in  their  radiation  time  intervals,  ~0.1  s  (due  to 
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5.  4.2.  Evolution  of  the  interaction  in  space  (altitude)  and  time.  Snapshots  of  the 
ctric  field,  optical  excitation  rates  of  the  5577  X  line  of  O  (lighter  line)  and  the  1st 
sitive  band  of  N2  (solid  line),  and  average  electron  energy  shown  for  the  80  to  180  /is 
le  interval  in  20  /is  increments  following  the  injection  of  a  Einj  =  20  V/m  EM  pulse 
70  km  under  ‘tenuous’  nighttime  conditions  (Figure  2.1  (a)). 
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quenching)  at  90  km  for  the  5577  X  line  and  ~6  /is  (due  to  spontaneous  radiation)  for 
the  1st  positive,  for  quenching  and  radiation  rates  see  Appendix  D. 

The  altitude  distribution  of  maximum  intensities  of  the  various  lines  stimulated  by  an 
EM  pulse  with  Einj  =  20  V/m  is  shown  in  Figure  4.3  (a).  Note  that  the  values  plotted 
here  represent  the  maximum  emission  intensity  reached  at  that  altitude  sometime  during 
the  260  /is  interaction  episode  which  was  illustrated  by  means  of  space-time  shapshots 
in  Figure  4.2.  The  emission  intensities  of  ‘allowed’  bands  integrated  over  the  100  /is 
duration  of  the  EM  pulse*  are  shown  in  Figure  4.3  (b),  whereas  the  temporal  variations 
of  the  height  integrated  (70-100  km)  intensities  are  shown  in  Figure  4.3  (c). 

The  altitude  distribution  of  maximum  intensities  of  the  various  lines  stimulated  by  an 
EM  pulse  with  Einj  =  25  V/m  is  shown  in  Figure  4.4  (a).  Note  that  the  values  plotted 
here  represent  the  maximum  emission  intensity  reached  at  that  altitude  sometime  during 
the  260  /is  interaction  episode.  The  emission  intensities  integrated  over  the  duration  of 
the  emissions  (~  100  /is  for  the  short  lasting  band  emissions,  ~0.5  ms  at  90  km  for  6300 
X  line  of  O,  and  ~0.1  s  at  90  km  for  5577  X  line  of  O)  are  shown  in  Figure  4.4  (b), 
whereas  the  temporal  variations  of  the  height  integrated  (70-100  km)  intensities  are  shown 
in  Figure  4.4  (c). 

For  an  EM  pulse  with  E{nj  =  10  V/m,  the  1st  and  2nd  positive  bands  of  N2  have  10 
to  12  times  weaker  maximum  and  time  integrated  intensities  than  those  for  a  pulse  with 
E,nj  =  25  V/m.  Emissions  from  molecular  ions  have  ~  1,000  times  weaker  maximum 
and  time  integrated  intensities.  Intensities  of  6300  and  5577  X  lines  of  O  are  lower  by 
a  factor  of  ~10  in  comparison  with  the  E{nj  =  25  V/m  case. 

Figure  4.5  shows  the  dependence  of  the  maximum  emission  rates  versus  the  electric 


*  Note  that  the  duration  of  the  optical  emissions  of  the  molecular  bands  is  also  ~100  /is,  although  the 
6300  and  5577  A  of  O  last  longer  ( — 0.5  ms  and  ~0.1  s  at  90  km  respectively). 
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Fig.  4.3.  Intensities  of  optical  emissions.  Top 
left  (a):  The  altitude  distribution  of  the  maxi¬ 
mum  intensities  produced  by  a  single  EM  pulse 
with  E{nj=  20  V/m  under  ‘tenuous’  nighttime 
conditions  (Figure  2.1  (a)).  Top  right  (b):  The 
altitude  distribution  of  emission  intensities  inte¬ 
grated  over  the  100  ps  duration  of  the  EM  pulse 
for  the  same  conditions.  Bottom  left  (c):  The 
temporal  variation  of  the  height  integrated  in 
tensities  for  the  same  conditions.  Dashed  lines 
show  night  airglow  background  levels  [ Cham¬ 
berlain, ,  1978;  pp.  214-218],  a  for  the  1st  pos¬ 
itive  band  of  N2  and  the  Meinel  band  of  NJ,  /? 
for  the  1st  negative  band  of  OJ,  the  6300  and 
5577  A  lines  of  O,  and  7  for  the  2nd  positive 
band  of  N2,  and  the  1st  negative  band  of  NJ. 


field  amplitude  of  the  pulse,  expressed  in  terms  of  £100*.  Also  indicated  are  the  percent¬ 
ages  of  occurence  of  CG  discharges  with  ^ioo  larger  than  the  corresponding  numbers  in 
V/m,  based  on  ground-based  measurements  [Krider  and  Guo ,  1983]. 

We  see  from  Figure  4.5  that  due  to  the  somewhat  different  dependence  of  the  emission 
intensity  of  the  various  lines  on  £100,  the  spectral  content  of  the  optical  emissions  from 
the  D  region  vary  in  predictable  manner.  In  particular,  we  note  that  the  1st  and  the  2nd 
positive  bands  of  N2  and  6300  and  5577  A  lines  of  O  do  not  increase  as  rapidly  with 

*  Note  from  Section  2.5  that  E\oo  =  E{n j 
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iq8  Fig.  4.4.  Intensities  of  optical  emissions.  Top 
left  (a):  The  altitude  distribution  of  the  maxi- 
mum  intensities  produced  by  a  single  EM  pulse 
10  with  Einj=  25  V/m  under  ‘tenuous’  nighttime 
conditions  (Figure  2.1  (a)).  Top  right  (b):  The 
10  4  altitude  distribution  of  emission  intensities  inte¬ 
grated  over  the  duration  of  the  emissions  (see 
102  text)  for  the  same  conditions.  Bottom  left  (c): 

The  temporal  variation  of  the  height  integrated 
1Qo  intensities  for  the  same  conditions.  Dashed 
lines  show  airglow  background  levels  at  night 
[ Chamberlain ,  1978;  pp.  214-218],  a  for  the 
1st  positive  band  of  N2  and  Meinel  band  of  NJ, 
/ 3  for  the  1st  negative  band  of  OJ,  the  6300  and 
5577  A  lines  of  O,  and  7  for  the  2nd  positive 
band  of  N2,  and  the  1st  negative  band  of  N£. 


increasing  E 100  as  the  emission  bands  of  the  ionized  molecules.  For  example,  the  ratio 
of  emission  rates  of  the  1st  positive  of  N2  and  the  1st  negative  of  N£  for  discharges  with 
£100  =  20  V/m  (stronger  pulses  occur  ~5  %  of  the  time)  and  E^oo  =  10  V/m  (stronger 
pulses  occur  ~40  %  of  the  time)  is  ~7  and  ~400  respectively.  This  circumstance 
results  from  the  generally  higher  excitation  thresholds  of  the  emissions  for  the  ions 
and  shows  the  potential  use  of  spectroscopic  measurements  of  airglow  enhancements 
as  a  means  to  possibly  assess  the  strength  of  the  EM  pulse,  or  if  that  is  well  known 
via  other  measurements,  other  quantities  such  as  the  ambient  electron  density  (also  see 
Figure  4.7  below). 
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Fig.  4.5.  Maximum  emission  rates  as  a  function  of  the  electric  field  amplitude  E ioo- 
The  calculations  were  performed  under  ‘tenuous’  nighttime  conditions  (Figure  2.1  (a)). 
Also  indicated  are  the  percentages  of  occurrence  of  CG  discharges  with  E joo  larger  than 
the  corresponding  numbers  in  V/m  [Krider  and  Guo ,  1983]. 


Although  not  shown  here,  altitude  and  time  integrated  intensities  show  similar  depen¬ 
dences  as  a  function  of  the  input  electric  field  amplitude  since  the  peak  intensities  tend 
to  dominate  the  averages  over  time  and  altitude. 


4.2  RESULTS  FOR  A  ‘DENSE’  NIGHTTIME  D  REGION 

For  a  denser  D  region  the  reflection  altitude  of  the  EM  pulse  is  lower  (see  Section  3. 1  for 
a  short  discussion  on  reflection  height),  resulting  in  higher  local  neutral  density  and  less 
electron  heating  compared  to  a  ‘tenuous’  D  region.  Thus,  generally  larger  electric  field 
amplitudes  are  needed  to  produce  the  same  levels  of  excitation  of  optical  emissions. 

As  a  ‘dense’  ambient  density  profile,  we  consider  a  profile  (b)  in  Figure  2.1,  which 
may  represent  nighttime  conditions  during  enhanced  geomagnetic  activity  or  in  the  auroral 
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regions.  In  accordance  with  the  above  discussion  the  amplitudes  of  EM  pulses  that  we 
consider  in  this  section  are  larger  than  in  Section  4.1  for  a  ‘tenuous’  D  region.  For  an  EM 
pulse  with  initial  intensity  Ein}  =  35  V/m  the  altitude  distribution  of  maximum  intensities 
is  shown  in  Figure  4.6  (a).  The  emission  intensities  integrated  over  the  duration  of  the 
emissions  (~100  /j, s  for  the  short  lasting  band  emissions,  ~0.5  ms  at  90  km  for  6300  A 
line  of  O,  and  ~0.1  s  at  90  km  for  5577  A  line  of  O)  are  shown  in  Figure  4.6  (b),  whereas 
the  temporal  variation  of  the  height  integrated  intensities  is  given  in  Figure  4.6  (c). 
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Fig.  4.6.  Intensities  of  optical  emissions.  Top 
left  (a):  The  altitude  distribution  of  the  maxi¬ 
mum  intensities  produced  by  a  single  EM  pulse 
with  Etnj  =35  V/m  under  ‘dense’  nighttime  con¬ 
ditions  (Figure  2.1  (b)).  Top  right  ( b ):  The 
altitude  distribution  of  emission  intensities  inte¬ 
grated  over  the  duration  of  the  emissions  for  the 
same  conditions.  Bottom  left  ( c ):  The  temporal 
variation  of  the  height  integrated  intensities  for 
the  same  conditions.  Dashed  lines  show  airglow 
background  levels  at  night  [ Chamberlain ,  1978; 
pp.  214-218],  a  for  the  1st  positive  band  of  N2 
and  the  Meinel  band  of  N£,  /?  for  the  1st  nega¬ 
tive  band  of  O2,  the  6300  and  5577  A  lines  of 
O,  and  7  for  the  2nd  positive  band  of  N2,  and 
the  1st  negative  band  of  N£. 
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For  ExnJ  =  30  V/m,  the  intensities  of  the  1st  and  2nd  positive  band  of  N2  are  ~1.5 
times  smaller  than  for  a  Einj  =  35  V/m  pulse.  Emissions  from  molecular  ions  have  2 
to  3  times  weaker  maximum  and  time  integrated  intensities.  Intensities  of  the  6300  and 
5577  A  lines  of  O  are  lower  by  a  factor  of  ~1.5  in  comparison  with  the  EinJ  =  35  V/m 
case. 

For  Einj  =  20  V/m,  the  intensities  of  the  1st  and  2nd  positive  bands  of  N2  have 
~10  times  weaker  maximum  and  time  integrated  intensities  than  those  for  a  Exnj  =  35 
V/m  pulse.  Emissions  from  molecular  ions  have  ~200  times  weaker  maximum  and  time 
integrated  intensities.  Intensities  of  the  6300  and  5577  A  lines  of  O  are  lower  by  a  factor 
of  ~10  in  comparison  with  the  E{nj  =  35  V/m  case. 

Figure  4.7  shows  the  dependence  of  maximum  emission  rates  on  the  electric  field 
amplitude  of  the  pulse,  expressed  in  terms  of  JE/joo- 

Also  indicated  are  the  percentages  of  occurence  of  CG  discharges  with  £100  larger  than 
corresponding  numbers  in  V/m,  based  on  ground-based  measurements  [Krider  and  Gun, 
1983].  The  dependences  on  JE'ioo  are  similar  to  the  ‘tenuous’  case  but  with  important 
quantitative  differences  in  the  the  relative  amplitudes  of  the  emission  lines  for  a  given 
.Eioo-  For  example,  for  E'ioo  =  20  V/m,  the  ratio  of  the  intensity  of  the  1st  positive  band  of 
N2  to  that  of  the  1st  negative  of  is  ~103  for  the  ‘tenuous’  case  (Figure  4.5)  but  r 2x10s 
for  the  ‘dense’  D  region  case  (Figure  4.7).  Comparison  of  Figures  4.5  and  4.7  show 
the  potential  application  of  spectroscopic  measurements  of  lightning-induced  enhanced 
airglow  to  determine  the  ambient  D  region  electron  density  profile. 

4.3  COMPARISON  WITH  EARLIER  WORK  [TARANENKO  ET  AL„  1992] 

In  this  section,  we  compare  our  results  with  those  presented  in  an  earlier  paper,  which 
was  based  on  a  simplified  model  of  the  interaction  neglecting  reflection  of  the  pulse 
and  assuming  a  Maxwellian  distribution  of  electrons  and  which  used  somewhat  different 
ambient  density  profiles. 
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Fig.  4.7.  Maximum  emission  rates  as  a  function  of  the  electric  field  amplitude  E\oo- 
The  calculations  were  performed  under  ‘dense’  nighttime  conditions  (Figure  2.1  (b)). 
Also  indicated  are  the  percentages  of  occurence  of  CG  discharges  with  E\oo  larger  than 
corresponding  numbers  in  V/m  [Krider  and  Guo ,  1983]. 


Results  presented  in  Figures  3  to  6  of  Taranenko  et  al.  [1992]  for  the  so-called  Case  A 
correspond  to  a  E\qq  =  20  V/m  pulse,  which  in  turn  corresponds  in  our  case  to  a  Einj  =  25 
V/m  pulse  (Section  2.5).  To  make  a  comparison  of  the  two  models  we  used  our  new 
model  described  in  this  thesis  for  Etnj  =  25  V/m  with  the  electron  density  profile  used  by 
[Taranenko  et  al.,  1992].  Results  indicate  that  the  kinetic  calculations  give  substantially 
lower  emission  intensities  for  the  same  EM  pulse  amplitude  compared  with  the  simple 
model  of  Taranenko  et  al.  [1992].  Namely,  the  maximum  intensity  of  the  1st  positive 
band  of  N2  is  5.2  x  107  cm-3  sec-1,  i.e.  30  times  lower.  The  corresponding  numbers 
for  other  lines  are:  the  2nd  positive  band  of  N2  is  90  times  lower,  Meinel  band  of  N£ 
is  500  times  lower,  the  1st  negative  band  of  N£  is  2,700  times  lower,  the  1st  negative 
band  of  O2  is  600  times  lower,  the  6300  A  line  of  O  is  25  times  lower,  and  the  5577  A 
line  of  O  is  30  times  lower.  These  differences  come  about  due  to  several  factors.  Firstly, 
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due  to  reflection  of  the  pulse,  the  maximum  heating  occurs  at  ~88  km  (5  km  lower  than 
in  the  old  model)  where  the  collision  frequency  is  much  higher,  resulting  in  less  heating 
(6  eV  average  energy  instead  of  9  eV)  and  therefore  a  reduced  number  of  electrons 
above  the  excitation  thresholds  of  most  emissions.  Secondly,  the  precipitous  drop  (due  to 
high  energy  losses  for  energies  above  the  thresholds  of  vibrational,  optical,  dissociative, 
and  ionizational  excitations;  see  Figure  2.3  (b))  of  the  electron  distribution  function 
decreases  the  number  of  energetic  electrons  in  the  tail  of  the  distribution  compared  with 
a  Maxwellian  distribution  of  the  same  average  energy. 

For  optical  emissions  with  low  excitation  thresholds  the  latter  factor  is  not  as  important. 
For  example,  for  the  5577  A  line  of  O,  the  5  km  decrease  in  the  altitude  of  maximum 
heating  results  in  ~5.75  times  less  electrons  and  ~2.75  times  less  oxygen  atoms,  and 
the  decrease  in  average  energy  gives  a  factor  of  ~2  decrease  in  intensity.  Alltogether 
this  results  in  ~31.6  times  lower  emission  rate  for  this  line  which  matches  the  30  times 
difference  between  the  new  and  old  models. 

We  also  note  also  that  in  [Taranenko  et  al.,  1992]  a  slightly  different  excitation  cross 
section  was  used  for  the  2nd  positive  band  of  N2  (see  Appendix  B).  The  usage  of  the 
different  cross  sections  does  not  play  an  important  factor  in  the  difference  between  the 
emission  intensities  reported  here  and  in  [ Taranenko  et  al.,  1992]. 


5 

Discussion 


In  this  chapter,  we  discuss  the  various  implications  of  our  results  and  compare  them  with 
available  experimental  data. 

5.1  ELECTRON  DENSITY  CHANGES  AND  ASSOCIATED  EFFECTS 


5.1.1  Perturbation  of  VLF  Subionospheric  Signals 

As  was  described  in  Chapter  1,  the  first  experimental  evidence  for  direct  upward  cou¬ 
pling  of  lightning  energy  to  the  lower  ionosphere  was  in  the  form  of  ‘early’  and/or  ‘fast’ 
perturbations  (sometime  called  ‘early’  or  ‘fast’  Trimpi  events)  of  subionospheric  VLF 
transmitter  signals  [Inan  et  al.,  1988;  1993].  Detailed  theoretical  investigation  of  such 
effects  needs  to  be  based  on  an  adequate  three  dimensional  model  of  the  interaction  to 
determine  three  dimensional  parameters  (most  importantly  the  transverse  extent)  of  the 
region  of  electron  density  changes  which  can  in  turn  be  used  in  a  computer-based  model 
of  subionospheric  VLF  propagation  [Poulsen  et  al.,  1993]  to  calculate  the  resultant  VLF 
amplitude  and  phase  changes.  Such  modeling  is  beyond  the  scope  of  our  present  one 
dimensional  model.  However,  some  first  order  estimates  of  the  expected  VLF  perturba¬ 
tions  due  to  the  electron  density  changes  produced  by  EM  lightning  pulses  as  calculated 
with  the  one  dimensional  model  can  be  made  on  the  basis  of  geometrical  factors. 
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Since  the  source  of  the  EM  waves  (i.e.,  lightning)  is  typically  located  at  a  distance 
of  about  80  to  90  kilometers  from  the  region  of  heating,  it  is  reasonable  to  assume  that 
its  radiation  pattern  is  not  well  directed,  hence  the  transverse  extent  of  the  distributed 
region  can  be  expected  to  be  of  the  same  order  as  the  distance  from  the  source  (80  to  90 
km).  On  this  basis,  [Inan  et  al.,  1991]  used  a  computer-based  model  of  subionospheric 
VLF  propagation  [Poulsen  et  al.,  1990]  to  calculate  the  amplitude  and  phase  changes 
due  to  electron  density  increases  produced  by  an  EM  pulse  ionizing  the  D  region.  They 
found  that  a  150  km  radius  ionization  patch  with  ~10  %  maximum  change  in  electron 
density  at  ~95  km  altitude  produces  amplitude  and  phase  changes  of  subionospheric  VLF 
waves  within  the  range  of  experimentally  observed  values  [Inan  et  al.,  19886].  Also, 
two  dimensional  modeling  of  heating  based  on  a  simple  model  assuming  a  Maxwellian 
electron  distribution  [Rodriguez  et  al.,  1992]  showed  widths  at  half-maximum  of  the 
ionized  regions  of  90  km  (for  a  E\oo  =  20  V/m  pulse  from  a  vertical  radiator  at  0  km 
altitude)  to  260  km  (for  E\oo  =  40  V/m),  consistent  with  the  arguments  above. 

In  comparing  our  results  with  those  of  Inan  et  al.  [1991],  we  see  that  although  the 
general  conclusions,  such  as  the  fact  that  lightning-induced  heating  of  the  electrons  leads 
to  impact  ionization,  are  borne  out,  there  are  striking  differences.  The  results  of  [Inan  et 
al.,  1991]  were  based  on  a  simplified  model  of  coupling,  which  assumed  a  Maxwellian 
distribution  to  be  maintained  during  the  heating  and  which  neglected  the  reflection  of 
the  EM  pulse  from  the  lower  ionosphere.  The  estimation  of  the  resulting  ionization  was 
also  handled  in  terms  of  a  first  approximation  by  assuming  that  secondary  production 
occurred  only  for  electrons  with  average  energy  >35  eV.  It  is  interesting  to  note  that 
the  relatively  simple  model  of  Inan  et  al.  [1991]  appears  to  have  roughly  predicted 
the  magnitude  of  the  density  enhancements  (a  few  percent  to  a  few  tens  of  percent) 
that  we  find  for  the  higher  altitudes.  However,  the  earlier  model  predicted  only  density 
enhancements  while  our  results  also  show  a  reduction  in  density  at  the  lower  altitudes 
due  to  attachment,  which  produces  a  general  sharpening  of  the  density  profile.  Also, 
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the  altitude  of  maximum  perturbation  is  found  to  be  somewhat  lower  than  in  the  earlier 
model  due  to  the  reflection  of  the  pulse  from  the  ionosphere. 

In  terms  of  the  relevance  of  our  results  to  the  interpretation  of  experimental  data,  we 
note  that  disturbances  at  lower  altitudes  and  with  a  more  drastic  change  from  ambient 
(i.e.,  both  an  increase  at  highter  altutudes  and  a  reduction  at  lower  heights)  would  in 
general  be  more  likely  to  lead  to  subionospheric  VLF  signal  perturbations.  Thus,  the 
predicted  modifications  in  electron  density  of  1-30  %  of  the  ambient  produced  by  single 
strokes  should  be  quite  sufficient  to  produce  detectable  amplitude  and  phase  changes  of 
subionospheric  VLF  signals.  In  this  sense,  our  results  reaffirm  the  conclusions  of  [Inan 
et  al.,  1991],  in  that  this  heating  and  ionization  mechanism  appears  to  explain  all  known 
features  of  the  ‘early’  and/or*  fast’  subionospheric  VLF  signal  changes  which  occur  at 
night  nearly  simultaneously  with  lightning  discharges. 

Our  results  show  that  the  phenomenon  of  lightning-induced  ionization  and  attachment 
occurs  predominantly  under  conditions  of  a  ‘tenuous’  D  region  which  are  realized  during 
nighttime.  Physically,  this  effect  is  due  to  the  fact  that  for  a  ‘dense’  D  region  the 
reflection  altitude  of  the  lightning  EM  pulses  is  located  at  lower  altitudes  where  the  mean 
free  path  of  electrons  is  short  so  that  more  of  the  energy  goes  into  excitation  of  the 
molecular  energy  levels  which  require  smaller  amounts  of  energy  than  ionization  and 
dissociative  attachment.  Due  to  the  exponential  dependence  of  the  electron  mean  free 
path  versus  altitude,  the  critical  amplitude  of  the  electric  field  above  which  appreciable 
perturbations  are  produced  is  a  strong  function  of  the  D  region  electron  density  profile. 
For  example,  in  accordance  with  the  results  of  Chapter  3,  during  daytime  (profile  (c) 
of  Figure  2.1)  even  the  strongest  individual  pulses  (initial  intensity  E{nj  =  50  V/m)  do 
not  produce  substantial  changes  in  the  electron  density  ( | <5 A^e / A^e  [  <  3  x  10-3),  whereas 
during  nighttime  (profile  (a)  of  Figure  2.1)  an  EM  pulse  with  Etn}  =  10  V/m  (JSjoo  —  8 
V/m)  is  expected  to  produce  measurable  electron  density  changes  of  |/>7Ve/./Ve|  — 

The  strong  dependence  of  the  interaction  on  ambient  electron  density  may  explain  recent 
nighttime  observations  of  ‘fast’/'early’  VLF  perturbation  events  [Inan  et  al.,  1993]  as 
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associated  with  CG  flashes  with  electric  field  amplitudes  as  low  as  E\qo  =  6  V/m,  as 
possibly  being  a  consequence  of  an  unusually  tenuous  D  region. 


5.1.2  Formation  of  Ionization  ‘Bubbles’  over  Thunderstorms 

An  important  prediction  of  [Inan  et  al.,  1991]  was  the  possible  formation  of  ionization 
‘bubbles’  over  thunderstorm  centers  due  to  the  accumulation  of  ionization  produced  by 
successive  strokes.  This  is  reminiscent  of  sporadic  E  layers  occasionally  observed  in 
the  lower  ionosphere  [e.g.,  Fenn  et  al.,  1987;  p.  21-55],  sometimes  in  association  with 
lightning  [Isted,  1954].  Our  kinetic  model  shows  that  this  is  indeed  the  case,  with  the 
added  caveat  that  the  density  is  depleted  at  lower  altitudes  and  that  the  density  profile 
is  substantially  sharpened  as  the  ‘bubble’  builds  up,  e.g.,  the  scale  height  decreases  by 
a  factor  of  ~2  at  ~85  km  for  a  single  Einj  =  20  V/m  EM  pulse.  An  example  of  such 
a  development  of  an  electron  density  disturbance  is  shown  in  Figure  5.1.  One  would 
expect  that  such  sharpened  electron  density  profiles  over  thunderstorms  would  increase 
the  reflectivity  of  VLF  waves  from  the  ionosphere  and  improve  VLF  transmission  along 
the  Earth-ionosphere  waveguide.  In  this  connection,  it  is  interesting  to  note  that  all 
reported  cases  of  ‘early’  VLF  perturbations  have  been  positive  amplitude  changes  [Inan 
et  al.,  1993],  as  expected  from  such  a  modification  of  the  density  profile  as  shown  in 
Figure  5.1. 

Increases  in  the  electron  density  over  a  thunderstorm  center  produced  in  the  initial 
stages  of  a  thunderstorm  may  possibly  lead  to  a  denser  D  region  and  prevent  production 
of  appreciable  changes  in  electron  density  by  the  lightning  EM  pulses  occuring  in  later 
stages.  Such  an  effect  may  occur  due  to  the  decrease  in  the  reflection  height  of  the 
EM  pulses  to  the  altitudes  where  most  of  their  energy  transfers  to  excitations  other  than 
dissociative  attachment  and  ionization.  Hence,  observed  characteristics  of  ‘early’  and/or 
‘fast’  VLF  perturbation  events  may  evolve  during  a  thunderstorm;  one  such  example  was 
observed  in  the  case  of  a  thunderstorm  center  which  stayed  active  and  centered  on  a  VLF 
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Electron  Density  (cm-3) 

Fig.  5.1.  Evolution  of  an  electron  density  disturbance  due  to  successive  EM  pulses 
from  lightning.  The  disturbances  were  produced  by  an  EM  wave  with  E,nj  =  25  V/m  and 
different  time  durations  of  0.1  ms  (a  single  pulse),  1  ms  (equivalent  to  a  sequence  of  about 
8  pulses),  1.6  ms  (equivalent  to  a  sequence  of  about  14  pulses),  and  2.2  ms  (equivalent 
to  a  sequence  of  about  20  pulses)  under  ‘tenuous’  nighttime  conditions  (Figure  2.1  (a)). 

path  for  ~12  hours  although  ‘early’  VLF  events  were  only  observed  during  the  first  two 
hours  [Inan  et  al.,  19886]. 

5.1.3  Stimulation  of  ULF  waves 

The  conductivity  of  a  collisional  magnetized  plasma  in  a  quasi-constant  electric  field  is 
an  anisotropic  function  of  the  electron  and  ion  gyrofrequencies  and  the  electron  collision 
frequency  [e.g.,  Rees,  1989;  p.  202].  Since  all  three  components  of  the  ionospheric 
conductivity  tensor  are  proportional  to  the  electron  density  Ne,  changes  in  the  ionospheric 
electron  density  can  modify  the  conductivity  and  lead  to  variations  in  ionospheric  currents 
[e.g..  Bell,  1976].  In  the  context  of  our  results,  slow  variations  of  the  electon  density 
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above  thunderstorms  would  be  expected  since  typical  flash  rates  [ Uman ,  1987,  p.  49]  are 
on  the  order  of  several  flashes  every  ten  seconds,  and  the  density  relaxation  time  is  on 
the  order  of  tens  of  seconds  [e.g.,  Glukhov  et  al.,  1992].  Such  slow  variations  occuring 
over  few  tens  of  second  time  scales  can  be  expected  to  produce  ionospheric  current 
perturbations  in  the  ULF  range.  Given  the  expected  large  horizontal  scale  (hundreds  of 
kilometers,  see  discussion  above  and  Rodriguez  et  al.  [1992])  of  the  disturbed  ionospheric 
regions,  these  current  perturbations  can  result  in  increased  ULF  wave  generation  in  the 
vicinity  of  thunderstorms.  Recent  experimental  results  [Fraser-Smith,  1993]  showing 
increased  ULF  activity  near  thunderstorms  may  be  experimental  evidence  for  such  effects. 


5.2  EXCITATION  OF  OPTICAL  EMISSIONS 

In  this  section  we  discuss  several  issues  related  to  remote  detection  and  measurement  of 
the  optical  signatures  of  lightning-induced  heating.  These  are:  (i)  emission  intensities 
compared  to  ambient  airglow  background;  and  (ii)  radiative  transfer.  We  also  compare 
our  model  results  directly  with  the  observations  of  Boeck  et  al.  [1992]. 

5.2.1  Emission  Intensities  Compared  to  Airglow  Background 

The  comparison  of  the  predicted  optical  emission  intensities  with  the  ambient  airglow 
background  is  critically  important  for  assessing  their  detectability  in  different  experi¬ 
mental  configurations.  The  intensities  of  emissions  for  vertical  observations  are  shown 
in  Figures  4.4  (c)  and  Figures  4.6  (c)  together  with  the  known  ambient  airglow  levels 
for  each  line  or  band  separately  indicated.  The  comparison  of  the  EM  pulse  stimulated 
emissions  with  the  ambient  levels  can  be  summarized  as  follows. 

The  intensities  of  the  6300  and  5577  A  lines  of  O  have  low  absolute  values  for  all 
cases  considered  here  and  are  more  than  one  order  of  magnitude  less  than  the  ambient 
intensity  of  the  night  airglow  for  zenith  observations,  which  is  ~102  R  [ Chamberlain , 
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1978;  p.  215]  for  vertical  observations.  This  fact  makes  these  emissions  difficult  to 
detect  in  practice. 

The  intensities  of  the  transitions  of  N2  (the  1st  and  2nd  positive  bands),  N£  (the 
1st  negative  and  Meinel  bands),  and  OJ  (the  1st  negative  band)  are  several  orders  of 
magnitude  larger  than  typical  background  emissions  in  corresponding  spectral  regions 
(see  Figures  4.3  (c),  4.4  (c),  and  4.6  (c)  and  [Chamberlain,  1978;  pp.  214-218])  for 
vertical  observations.  Hence,  all  of  the  molecular  band  emissions  stimulated  by  EM  pulses 
from  lightning  discussed  in  this  thesis  are  substantially  above  the  background  levels  at 
night  and  can  be  registered  using  specially  designed  high  sensitivity  photometers. 

In  terms  of  detectability  with  a  typical  auroral  airglow  photometer,  we  note  that  a 
good  photometer  would  have  a  sensitivity  of  ~50  Rayleigh-seconds  (R-s)  [5.  Mende, 
private  communication].  In  this  context  we  note  that  the  1st  and  2nd  positive  bands  of 
N2  have  intensities  of  1,200  and  200  R-s  and  should  be  detectable  from  the  ground  with 
conventional  instruments.  However,  other  emissions  have  lower  than  50  R-s  intensities 
and  would  have  to  be  observed  with  specially  designed  instruments.  Such  considerations 
are  important  in  the  design  of  optical  instruments  aimed  at  detecting  lightning-induced 
emissions  from  the  D-region. 

The  fact  that  the  longest  radiation  time  of  the  strong  allowed  emissions  considered 
above  is  ~14  /xs  indicates  that  relatively  fast  optical  systems  with  50-100  /xs  (about  the 
time  interval  of  emissions  from  a  single  EM  lightning  pulse)  to  several  seconds  time 
of  exposure  should  be  used  in  experiments  aimed  at  observing  the  optical  emissions 
from  the  D-region  associated  with  lightning.  In  this  context,  we  note  that  the  intense 
radio  impulse  from  lightning  can  be  used  to  automatically  trigger  the  camera.  Since 
the  radiation  time  of  the  brightest  emissions  is  very  small,  it  is  important  to  take  into 
account  the  size  of  the  emitting  region  when  one  calculates  the  intensity  and  duration  of 
the  remotely  observed  optical  pulse.  This  consideration  is  particularly  important  for  limb 
observations  and  all-sky  vertical  observations  and  in  view  of  the  fact  that  the  size  of  the 
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observed  emitting  region  can  reach  500  km  [Boeck  et  al.,  1992],  In  such  cases,  the  light 
travel  delay  would  spread  out  the  duration  of  the  impulse  from  ~50  /is  to  about  1 .5  ms 
and  consequently  lower  its  intensity  as  compared  to  the  case  when  the  light  travel  time 
is  negligible  in  comparison  with  the  duration  of  the  emissions. 

5.2.2  Radiative  Transfer 

Radiative  transfer  of  optical  emissions  in  the  atmosphere  (i.e.,  the  process  of  photon 
propagation  in  the  atmosphere  and  processes  of  scattering,  absorption,  and  re-emission 
that  occur  during  such  propagation)  is  important  for  remote  sensing  of  the  heated  region, 
especially  from  ground-based  platforms.  The  rigorous  treatment  of  these  processes  is  an 
area  of  active  research  by  itself  (see  for  example  [Fenn  et  al,  1985;  Chapter  18]),  so  in 
this  section  we  merely  make  order  of  magnitude  estimates  that  show  optical  extinction 
in  the  atmosphere  for  each  band  or  line  considered. 

The  1st  and  2nd  positive  bands  of  N2  are  not  resonantly  absorbed  in  the  atmosphere. 
They  originate  from  transitions  between  excited  levels  (B3irg  — »  A3!^  and  C3ttu  — ►  B3tt9 
transitions  respectively,  as  shown  in  Figure  4.1)  with  a  minimum  energy  of  ~6  eV. 
Hence,  resonant  absorbers  for  these  emissions  are  practically  not  present  in  the  cool 
(~250°  K)  atmosphere.  In  principle,  these  emissions  could  be  resonantly  absorbed  in  the 
heated  region;  however,  the  region  is  optically  transparent  due  to  low  density  of  excited 
molecules  (even  in  the  extreme  case  of  up  to  109  excited  particles  per  cc,  with  a  cross 
section  of  ~  10-16  cm2,  and  a  106  cm  path  we  have  optical  depth  of  only  0.1,  i.e.  the 
radiation  is  attenuated  by  a  factor  of  only  exp(— 0. 1)  ~  0.9  after  passing  through  the 
region).  Rayleigh  scattering  (i.e.,  scattering  by  induced  dipole  moment  of  molecules)  is 
the  major  mechanism  that  determines  the  optical  thickness  of  the  clear  atmosphere  for 
these  bands  and  it  results  in  intensity  losses  of  up  to  4  %  for  the  red  and  up  to  25  %  for 
the  blue  spectral  ranges  per  vertical  traverse  through  the  atmosphere  [Fenn  et  al.,  1985]. 
Rayleigh  scattering  loss  decreases  to  a  few  percent  for  observations  from  a  high  altitude 

20  km)  aircraft. 


DISCUSSION 


67 


Meinel  and  1st  negative  bands  of  N£  and  the  1st  negative  band  of  O2  also  originate 
from  transitions  between  excited  levels  with  a  minimum  energy  of  ~15  eV  [Vallance 
Jones,  1974;  p.  90]  so  that  the  same  arguments  discussed  above  for  the  N2  bands  also 
apply  to  them. 

5.2.3  Comparison  with  Boeck  et  al.  [1992]  Observations 

We  now  show  that  recent  observations  of  lightning  induced  brightening  in  the  airglow 
layer  from  the  Space  Shuttle  [Boeck  et  al,  1992]  can  be  explained  as  electron  impact 
induced  emissions  from  the  heated  region  above  a  powerful  lightning  stroke,  produced 
in  the  manner  described  in  this  thesis. 

Long  lasting  emissions  from  O  and  O2  [Taranenko  et  al.,  1992]  were  not  detected  on 
the  Space  Shuttle  for  two  reasons:  (1)  the  brightest  of  them,  Herzberg  I  band  of  O2  lies 
outside  of  the  camera  wavelength  sensitivity  region,  (ii)  the  expected  excitation  levels 
of  the  6300  A  and  5577  A  lines  of  O,  and  the  atmospheric  band  of  O2  is  below  the 
sensitivity  of  the  camera  [ Taranenko  et  al.,  1992]. 

Each  of  the  short  lasting  bands  discussed  above  has  transitions  with  wavelengths  from 
3600  to  7200  A  to  which  the  TV  camera  on  the  Space  Shuttle  was  sensitive  with  a 
peak  sensitivity  at  4400  A.  However,  transitions  between  states  with  low  vibrational 
numbers  of  the  1st  and  2nd  positive  bands  of  N2,  and  Meinel  band  of  N2  have  their 
wavelengths  outside  of  the  sensitivity  region  of  the  camera.  Only  two  strong  lines  of 
the  1st  negative  band  of  N|  (0-0  transition  at  3914.4  A  and  0-1  transition  at  4278.1  A) 
lie  near  the  peak  of  the  camera  sensitivity.  To  obtain  the  intensity  of  the  emissions  for 
nighttime  limb  observations  in  Rayleighs  one  must  multiply  numbers  on  the  ordinate  scale 
of  Figure  4.4  by  the  thickness  of  the  emitting  region  given  in  units  of  10  km  (e.g.,  for  a 
500  km  region  as  given  in  [Boeck  et  al.,  1992]  the  multiplier  is  50).  As  is  evident  from 
Figure  4.4  (c),  the  time  duration  of  the  emissions  for  a  100  /xs  EM  pulse  with  Ein}  =  25 
V/m  is  ~50  /xs  (resulting  in  15  km  thickness  of  the  emitting  region  of  simultaneously 
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registered  photons)  so  that  the  side  view  intensity  of  this  band  averaged  over  5  km  of 
altitude  is  ~5x  104  R.  The  observed  duration  of  the  emissions  for  the  500  km  thickness 
of  the  emitting  region  is  ~  1 .5  ms.  Hence,  the  observed  effective  intensity  in  R-s  at  the 
Space  Shuttle  would  be  ~10  times  lower  due  to  1/60  s  integration  time  of  the  camera. 
After  this  reduction,  the  intensity  observed  at  the  Shuttle  would  be  ~5x  10^  R.  Emissions 
of  the  1st  negative  band  of  Oj  are  in  the  spectral  region  of  the  camera  sensitivity  (~5200- 
6500  A).  They  have  similar  intensity  and  time  duration  to  those  of  the  1st  negative  band 
of  Nj  and  make  a  comparable  contribution  to  the  total  effective  intensity  ~104  R  of  the 
airglow  transient  over  the  visible  spectrum  range.  For  the  E{nj  =  20  V/m  case  at  night 
the  effective  intensity  of  the  emissions  is  ~10  times  lower  than  for  the  case  Exnj  =  25 
V/m  (Section  4.3). 

As  has  been  mentioned  in  section  1 .3,  the  intensity  of  the  airglow  transient  observed 
from  the  Shuttle  [Boeck  et  al.,  1992]  was  about  twice  that  of  the  ambient,  which  is  several 
kR  over  the  visible  spectral  range  at  night  [Chamberlain,  1978;  pp.  214-218].  Hence,  in 
accordance  with  the  above  estimates  only  EM  pulses  with  Etnj  >  25  V/m  could  produce 
such  intensities.  Emissions  produced  by  EM  pulses  with  Einj  <  20  V/m  would  not  be 
distinguishable  from  the  background,  with  the  TV  camera  used  in  the  Boeck  et  al.  [1992] 
experiments. 

To  assess  the  lightning  intensity  that  would  lead  to  optical  emissions  at  the  detected 
level  [Boeck  et  al.,  1992],  we  note  that  E\oo  corresponds  to  Einj  x  (see  Sec¬ 

tion  2.5  for  discussion),  hence  E{nj  =  25  V/m  corresponds  to  jEioo  =  20  V/m.  According 
to  Krider  and  Guo  [1983]  CG  discharges  with  such  intensities  are  expected  to  occur 
~5  -  10  %  of  the  time.  Hence,  one  can  conclude  that  the  observed  airglow  transient 
[Boeck  et  al.,  1992]  was  due  predominantly  to  emissions  of  the  1st  negative  bands  of 
Nj  and  Oj  excited  by  a  relatively  intense  (and  thus  relatively  rare)  lightning  discharge. 
The  intensity  of  such  emissions  is  a  strong  function  of  the  EM  field  of  the  pulse  due  to 
the  high  excitation  thresholds  (~18.8  eV  and  ~18.2  eV  respectively)  and  is  expected  to 
be  below  the  background  level  at  night  for  the  EM  pulses  with  jEioq  <  16  V/m.  This 
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realization  is  consistent  with  the  fact  that  the  airglow  enhancement  reported  by  Boeck 
et  al.  [1992]  was  only  observed  once  during  the  period  of  several  hours  of  observation 
from  the  Space  Shuttle  (see  Section  1.2)  during  that  mission. 


6 

Summary  and  Suggestions  for  Future 
Research 


6.1  SUMMARY 

The  thrust  of  this  work  has  been  the  simulation  of  the  interaction  of  lightning  radiated  EM 
pulses  with  the  lower  ionosphere  by  means  of  a  self-consistent  solution  of  the  Boltzmann 
kinetic  equation  for  electrons  and  Maxwell’s  equations  for  the  electromagnetic  (EM) 
fields.  Our  results  indicate  that  attachment,  ionization,  and  optical  emissions  produced 
in  this  interaction  are  sufficient  to  explain  all  heretofore  known  features  of  the  ‘early’ 
and/or  ‘fast’  VLF  signal  perturbations  and  airglow  brightening  observed  simultaneously 
with  lightning.  Our  model  also  provides  the  basis  for  quantitative  evaluation  of  a  wide 
variety  of  other  electrodynamic  and  chemical  processes  (especially  those  that  deal  with 
the  altered  electron  and  negative  ion  densities)  in  the  ionosphere  which  might  occur  in 
response  to  the  EM  energy  released  in  the  lightning  discharges. 

Under  nighttime  conditions,  individual  pulses  of  10-20  V/m  (normalized  to  100  km 
distance)  produce  changes  in  electron  density  of  1-30  %  of  the  ambient  while  a  sequence 
of  pulses  leads  to  >100  %  modification  at  altitudes  <95  km.  The  nature  of  the  density 
changes  produced  by  lightning  radiation  are  such  that  they  lead  to  a  ‘sharpening’  of 
the  lower  ionospheric  boundary  by  causing  a  reduction  in  electron  density  at  75-85  km 
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(due  to  attachment)  and  a  substantial  increase  at  85-95  km  (due  to  ionization),  e.g.,  the 
scale  height  decreases  by  a  factor  of  ~2  at  ~85  km  for  a  single  20  V/m  EM  pulse  (see 
Figure  5.1).  The  slow  variations  in  the  electron  density  resulting  from  ionization  due 
to  succesive  flashes  should  lead  to  variations  in  the  electrical  conductivity  of  the  lower 
ionosphere  which  in  turn  might  also  stimulate  ULF  waves.  Also,  intense  electron  heating 
in  the  D  region  over  thunderstorms  causes  production  of  electron  density  ‘bubbles’  and 
may  significantly  alter  the  chemistry  in  that  region. 

Among  the  optical  emissions  the  most  promising  for  observations  are  determined  to 
be  the  1st  (red)  and  2nd  (blue)  positive  bands  of  N2,  which  emit  at  rates  of  ~7xl07 
and  ~107  cm'V1  at  ~92  km  respectively  for  a  £100  =  20  V/m  EM  pulse.  The  height 
integrated  intensities  of  the  emissions  reach  ~3x  107  and  ~6x  106  R,  lasting  for  ~50  /is 
and  producing  ~  1,200  and  ~200  Rayleigh-seconds  for  zenith  observations. 


6.2  SUGGESTIONS  FOR  FUTURE  RESEARCH 


6.2.1  Dependence  of  the  Lightning  Generated  Electron  Density  Changes  and 

Optical  Emissions  on  the  Temporal  Shape  and  Duration  of  Lightning  EM 
Pulses 

An  immediate  application  of  our  one  dimensional  model  is  to  extend  it  to  a  broader  range 
of  parameters  of  EM  pulses  from  lightning. 

As  a  first  step,  the  current  version  of  the  computer  code  can  be  applied  to  study  the 
response  of  the  ionosphere  to  a  variety  of  EM  pulse  shapes  for  CG  and  IC  discharges 
(see  Figure  1.1  for  examples  of  pulse  shapes).  As  a  result  of  such  a  study  some  spe¬ 
cial  features  of  the  pulse  shape  that  may  enhance  or  suppress  the  interaction  may  be 
revealed.  For  example,  the  constructive  interference  that  was  seen  to  be  very  effective  in 
Figure  3.1  would  not  occur  for  EM  pulses  which  were  dominantly  positive  or  negative 
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(rather  than  alternating),  thus  leading  to  smaller  perturbations  even  though  the  total  EM 
energy  of  the  pulse  may  be  the  same. 

Another  possible  way  to  extend  the  model  is  the  accomodation  of  larger  amplitudes 
of  injected  EM  pulses.  This  requires  some  modification  of  the  code  since,  due  to  the 
particular  numerical  scheme  employed,  the  present  version  has  an  upper  limit  for  the 
electric  field  amplitude  that  does  not  allow  the  simulation  of  interactions  with  the  lower 
ionosphere  of  EM  pulses  with  amplitude  Einj  >  25  V/m  for  a  ‘tenuous’  D  region. 
Although  EM  pulses  with  such  intensities  occur  less  than  5  %  of  the  time,  quantitative 
determination  of  the  extent  of  their  interaction  would  be  very  interesting,  especially  in 
view  of  the  highly  nonlinear  nature  of  coupling,  as  is  evident  from  Figure  3.4.  In  this 
context,  we  note  that  lightning  EM  pulses  with  intensities  of  as  much  as  i^ioo  ~100  V/m 
are  occasionally  observed  [e.g.,  Turman,  1977]. 

Investigation  of  the  interaction  of  very  short  (<10  fis)  EM  pulses  with  the  lower 
ionosphere  requires  the  removal  of  the  quasi  steady  state  assumption  from  the  model. 
Such  a  development  is  essentially  straightforward  but  would  require  two  to  three  orders 
of  magnitude  more  computer  memory  and  time  to  perform  the  computations  that  are 
necessary. 

6.2.2  Application  of  Model  Results  to  Data  on  Subionospheric  VLF  Signal 
Perturbations 

Thorough  theoretical  investigation  of  observed  ‘early’  VLF  perturbations  requires  an 
adequate  three  dimensional  model  of  the  interaction  to  determine  the  three  dimensional 
parameters  of  electron  density  changes  which  can  be  used  in  a  computer-based  model 
of  subionospheric  VLF  propagation  [Poulsen  et  al.,  1990]  to  calculate  the  amplitude  and 
phase  changes  due  to  the  electron  density  changes.  However,  some  first  order  estimates 
of  the  expected  VLF  perturbations  due  to  electron  density  changes  produced  by  lightning 
EM  pulses  can  be  made  on  the  basis  of  geometrical  factors  as  described  in  Section  5.1.1. 
Using  such  estimates,  calculations  similar  to  those  using  a  computer-based  model  of 
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subionospheric  VLF  propagation  [Poulsen  et  al.,  1993]  can  be  made  to  calculate  the 
amplitude  and  phase  changes  due  to  D  region  electron  density  changes  produced  by 
lightning  EM  pulses. 

6.23  Three  Dimensional  Modeling  of  the  Interaction  of  Lightning  EM  Pulses 
with  the  Ionosphere 

Another  obvious  future  development  of  our  one  dimensional  model  is  to  extend  it  to 
three  dimensions  using  the  quasi  steady  state  distribution  function  technique  for  three 
dimensional  values  of  electric  currents,  ionization  and  optical  emission  rates. 

We  expect  that  the  dimensions  of  the  disturbed  electron  density  region  will  differ  from 
that  of  the  enhanced  airglow  region  and  will  depend  on  the  discharge  strength,  discharge 
orientation,  dip  angle,  and  ionospheric  parameters.  Also,  the  3-D  modeling  would  allow 
us  to  assess  the  best  means  of  remote  sensing  of  the  regions  among  possible  techniques 
such  as:  VLF  remote  sensing  using  earth-ionosphere  waveguide  signals;  optical,  infrared, 
or  ultraviolet  remote  sensing  of  the  enhanced  airglow;  radio  interferometry  using  satellite 
beacons  or  natural  radio  sources;  or  ground-to  satellite  or  ground-to-rocket  radio  beacons. 

Results  of  the  three  dimensional  modeling  can  provide  detailed  inputs  to  facilitate 
quantitative  predictions  for  VLF  subionospheric  wave  perturbations  [Poulsen  et  al.,  1990], 
D  region  aeronomy  models  [Glukhov  et  al.,  1992],  ULF  wave  generation  [Bell,  1976], 
and  for  radiative  transfer  codes  [Fenn  et  al.,  1985,  p.  1 8-37]  that  would  describe  optical 
signatures  of  the  heated  region  as  observed  from  a  distance.  Such  modeling  can  also  allow 
the  determination  of  the  total  3-D  deposition  of  EM  energy  in  the  lower  ionosphere  and 
the  dominant  sinks  for  this  energy. 
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6.2.4  Possible  Effects  of  the  Lightning  EM  Pulses  on  the  Thermal  and  Chemical 
Balance  of  the  Upper  Atmosphere 

Our  results  concerning  the  interaction  of  individual  lightning  EM  pulses  with  the  lower 
ionosphere  bring  to  the  fore  the  question  of  the  significance  of  these  effects  on  a  global 
scale,  in  terms  of  the  integrated  effects  of  thunderstorms  worldwide  on  the  ionosphere. 
The  buildup  of  an  electon  density  ‘bubble’  in  the  lower  ionosphere,  expected  from  a 
sequence  of  EM  pulses  from  a  thunderstorm  were  mentioned  in  Section  5.1.2. 

In  further  studies,  the  Glukhov  et  al.  [1992]  model  of  the  D  region  chemistry  can  be 
applied  to  follow  the  evolution  of  such  an  electron  ‘bubble’.  Preliminary  applications  of 
such  a  model  to  some  previously  modeled  disturbances  indicate  that  some  of  the  least 
known  chemical  reaction  rates  in  the  nighttime  D  region  altitudes  may  be  measurable 
using  subionospheric  VLF  data  and  in  some  cases  new  ionization  layers  are  formed 
that  represent  a  different  attachment-detachment  quasi-equilibrium  value  from  that  of 
the  unperturbed  ambient  [Pasko  and  Inan,  1993].  Such  new  layers  may  exist  for  up  to 
~  105  s  following  the  disturbance.  This  is  reminiscent  of  sporadic  E  layers  occasionally 
observed  in  the  lower  ionosphere  [e.g.,  Fenn  et  al.,  1987;  p.  21-55]. 

There  are  some  other  indications  that  the  lightning  stimulated  processes  might  change 
the  lower  ionosphere  chemistry.  For  example,  a  three-body  collision  between  electrons 
and  two  oxygen  molecules 


e  +  02  +  02  ->  Oj  +  02  +  0.45eV  (6.1) 

is  considered  so  far  as  the  most  important  source  of  negative  ions  in  the  D  region  [Banks 
and  Kockarts,  1973].  However,  according  to  our  results  presented  in  Section  3.1  Fig¬ 
ure  3.3  (b),  EM  radiation  from  lightning  can  create  very  substantial  (on  the  order  of  50 
%  of  the  electron  density)  depletion  of  the  electron  density  due  to  dissociative  attachment 
of  electrons  to  O2 


e  +  02  — >  CT  +  O  -  7eV 


(6.2) 
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indicating  that  the  population  of  negative  atomic  oxygen  ions  O-  in  the  D  region  above 
a  thunderstorm  can  reach  up  to  50  %  of  the  positive  ion  density  in  that  region.  Such 
production  of  negative  ions  by  lightning-induced  heating  would  thus  effectively  compete 
with  reaction  (5.2)  in  producing  negative  ions. 

Once  created,  O-  is  transformed  into  other  negative  ions  by  means  of  charge  exchange, 
interchange,  clustering,  and  cluster  interchange  reactions,  involving  a  large  variety  of 
chemical  reactions  [Banks  and  Kockarts,  1973;  p.  281].  Hence,  appearance  of  negative 
atomic  oxygen  ions  O-  in  the  D  region  above  a  thunderstorm  can  be  expected  to  influence 
many  channels  of  chemical  reactions  in  that  region. 

6.2.5  Stimulation  of  ULF  waves 

As  was  discussed  in  Section  5.1.3,  changes  in  the  electron  density  in  the  lower  ionosphere 
modify  the  natural  conductivity  tensor  and  therefore  can  generate  ULF  waves.  To  evaluate 
the  importance  of  ULF  waves  generated  in  such  a  way  one  must  make  calculations 
of  ULF  currents  and  radiated  ULF  fields  stimulated  by  the  electron  density  changes 
and  make  comparison  of  the  radiated  fields  to  the  natural  magnetospheric  and  Earth- 
ionosphere  waveguide  backgrounds,  and  to  ULF  waves  stimulated  by  other  processes, 
such  as  electrons  precipitated  from  the  magnetosphere  by  whistlers  [Bell,  1976]. 


A 

Derivation  of  the  Truncated  Boltzmann 
Equations 


In  our  derivation  of  (2.4)  and  (2.5)  we  follow  the  general  pattern  outlined  in  Gurevich 
[1978,  p.  34].  We  start  with  general  form  of  the  Boltzmann  kinetic  equation  for  the 
electron  distribution  function  f(r,  v ,  t ) 


Of  _  df  e  r  ~  l 

~^7  +  v  ’  - 1  E  +  - 

at  or  my  c 


t7xB]}.|  +  5  =  ° 


(A.l) 


where  S  is  the  Boltzmann  collision  integral,  which  describes,  in  our  particular  case,  the 
change  of  the  electron  distribution  function  /  under  the  influence  of  electron  collisions 
with  neutrals.  The  general  expression  for  S  is  as  follows 


s  =  J  f  dvxdOq{u,  9)u{f(v)F(v\)  -  )}  (A.2) 

where  is  the  velocity  of  the  neutral  that  collides  with  the  electron,  u  =  |v  —  t>i|, 
q(u,0)  is  the  differential  effective  scattering  cross  section,  v 1  and  are  the  velocity  of 
the  electron  and  of  the  neutral  prior  to  their  collision  (their  velocity  after  the  impact  are 
respectively  v  and  uj);  and  F  is  the  distribution  function  of  the  neutrals.  The  integration 
in  equation  (A.2)  is  over  the  velocities  of  the  neutrals  (dv\)  and  over  the  scattering  angles 
dQ.  =  sin  9d0dp,  where  0  is  the  angle  between  the  vectors  (v  —  v\)  and  ( x?  —  tTj).  All 
other  variables  in  (A.l)  and  (A.2)  were  specified  in  Section  2.2. 
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The  first  simplification  which  we  made  in  writing  (A.  1)  deals  with  the  spatial  gradients 

a 

w)  °f  ^e  distribution  function;  we  neglect  such  gradients  since  they  are  only  important 
in  the  D  region  over  time  intervals  longer  than  several  hundred  seconds  [Tomko  et  al., 
1980]. 

If  the  electrons  are  under  the  influence  of  a  weak  electric  field  (the  specific  conditions 
are  given  below,  equation  (A.  11)),  such  that  their  random  velocity  is  much  larger  than 
their  average  directional  velocity,  then  the  distribution  function  depends  mainly  on  the 
absolute  value  of  the  velocity  and  not  on  its  direction.  Under  such  an  assumption,  the 
Boltzmann  equation  can  be  simplified  using  an  expansion  of  the  distribution  function  in 
spherical  functions  of  zero  order  (i.e.  Legendre  polynomials  Pk(cos0)  [Hochstrasser, 
1972])  in  velocity  space  [Allis,  1956]. 

For  the  sake  of  simplicity  in  the  following  analysis,  we  neglect  the  external  magnetic 
field.  A  similar  derivation  in  the  presence  of  a  constant  ambient  magnetic  field  was 
carried  out  by  [Davydov,  1937].  Accordingly,  we  consider  a  uniform  collisional  plasma 
in  the  presence  of  an  electric  field  E  parallel  to  the  2  axis  [Gurevich,  1978;  p.  36],  so  that 
there  is  only  one  preferred  direction  in  the  medium,  i.e.,  along  the  2  axis.  Consequently 

OO 

f(v,  r,  t)  =  ^  Pk(cos  9)fk(v ,  2,  t)  (A3) 

k=0 

where  0  is  the  angle  between  the  velocity  vector  and  the  2  axis.  Substitution  of  this 
expantion  into  equation  (A.l)  gives 


'£{P„(cos»)?A  -  ^cos«Pt(cos«)|£  -  +  5  =  0  (A.4) 


The  last  two  of  the  terms  inside  the  curly  brackets  result  from  the  fact  that 


e  d-l  =  E?i  -E^,edf  i  gsin^ 

dv  dvz  dv  v  d(cos  6) 
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Multiplying  equation  (A.2)  by  P*(cos0)  and  integrating  over  the  spherical  solid  angle 
dQ  =  sin  6d6d<p  we  obtain  the  following 


dh 

dt 


eEr  k  + 1  d  k+2f 
m  1(2/:  +  3)vk+2  dv  *+1 


+ 


kvk  1  d 
(2 k  -  l)3u 


+  Sk  =0 


(A. 5) 


Sk  =  —7—  [  SPk(cos0)dQ. 

4tt  J 

Where  we  have  taken  into  account  the  orthogonality  of  the  polynomials  [Hochstrasser, 
1972] 

J  ^  Pk(x)Pi(x)dx  =  2Y+\^k‘ 

(where  t>ki  =1  for  k  =  l  and  zero  otherwise)  and  the  recurrence  relations 


xPk 


k  + 1  k 

2k  +  1  Pk+]  +  2k  +  1 


(1  -x2) 


dl\ 

dx 


k(k  +  1) 

2  k  +  T 


[•Pjt-i  -  A+i) 


Equations  (A. 5)  represents  a  chain  of  equations  for  /* 


dfo 

dt 


e_E  d 

3mv2  dv 


(v2fi)  +  Sq  =  0 


(A. 6) 


dfi 

dt 


eE  \dfo  +  2  d 


m 


dv  5v 3  dv 


(v3h)\  +S\  =0 


dh 

dt 


eE  r2 


m 


rl&hUyi »]♦*-» 


etc.  This  chain  can  be  terminated  after  the  first  two  equations  if 


J2_d_ 

5v 3  dv 


(v3h) 


(A. 7) 

(A. 8) 


(A. 9) 
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We  neglect  fo  in  (A. 8)  and  using  approximations* 


dfi 

dt 


=  -*w/i, 


dfi 

dt 


=  -iuf2 


S  i  ^  vf\ ,  S2  ot  v/2 


where  i/  is  the  collision  frequency  determined  at  the  average  electron  energy  e.  We  obtain 
from  equations  (A. 8)  and  (A. 9) 


/i  = 


eE 


dfo 


m(—iu  +  v)  dv  ’ 


/2  = 


2eEv 


+ 1/) 


eE  3/o 


m(—iu>  +  v)v  dv 


The  condition  for  terminating  the  chain  is  therefore 


\df  >d  \  I  4  ^  ^  r  eE  d/ol'i. 

0  ' \5v^  dv{  m(—iu>  +  v)  dv  +  v)v  dv  .  ) 


(A.  10) 


For  the  average  electron  velocity  v,  we  have  w<i/|  and  assuming  ^  ~  l/t4  (A.  10) 
reduces  to 


e2 E2 /mev2  ■<  1 


(A.  11) 


Note  that  this  condition  is  satisfied  for  the  problem  at  hand  and  for  the  parameters  of 
interest.  For  example,  at  90  km  and  for  E  =  10  V/m,  e  =  4. 1  eV  (see  Figure  2.4), 
v{e)  =  107  s-1  so  that  the  left  hand  side  of  (A.  11)  is  ~  0.05  <C  1. 


In  the  presence  of  a  constant  magnetic  field,  the  same  procedure  can  be  applied  with 
the  additional  term  accounting  for  the  presence  of  the  magnetic  field  [Davydov,  1937] 


dfp  _  e  d  ( 
dt  3mv2  dv 


(v2E  ■  /l)  [v26v„ (1^  +„/„)]  +  S„,„  (A. 12) 


T  dfo 


*  These  approximations  are  equivalent  to  an  assumption  that  the  higher  order  expansion  terms  of  the 
distribution  function  vary  with  the  same  time  scale  (~l/u;)  as  the  electric  field. 

t  Holds  true  for  EM  pulses  from  lightning  with  VLF  type  spectra  in  the  lower  ionosphere  (see  Ap¬ 
pendix  C  for  numerical  values  and  additional  discussion). 

t  This  approximation  holds  true  for  a  Maxwellian  distribution  for  v  ~  tJ,  and  is  used  only  for  illustrative 
purposes  here. 
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df]  _  cE  dfo  e 
di  m  dv  me 

The  influence  of  the  external  magnetic  field  on  the  motion  of  energetic  electrons  in  the 
lower  ionosphere  is  weak  (see  Section  3.4  for  details)  and  the  equations  (A.  12)  and  (A.  13) 
are  valid  to  first  order  under  the  same  conditions  as  (A.  10)  or  (A.  11). 


B  x 


/.  -  v7\ 


(A.  13) 


B 

Electron  Impact  Excitation  Cross  Sections 


Plots  of  the  electron  impact  cross  sections  for  the  excitation  of  the  chosen  optical  emis¬ 
sions  as  a  function  of  electron  energy  are  shown  in  Figure  B.l.  The  data  on  electron 
impact  excitation  from  the  ground  states  compiled  in  this  figure  were  taken  from  the 
following  sources:  for  ]D  and  'S  states  of  O  [Doering  and  Gulcicek ,  1989],  for  i'E* 
and  A3E*  states  of  O2  [Rees,  1989,  Appendix  4  and  references  therein],  for  C37ru  and 
2?37rff  states  of  N2  [Cartwright  et  al.  1977]  with  30  %  reduction  in  accordance  with 
Phelps  [1987]  for  the  C'37ru  state,  for  A2nu  and  B2I.+  states  of  N£  [Cartwright  et  al. 
1975],  for  &4I"  state  of  0£  [Watson  et  al.  1967], 

Plots  of  the  electron  impact  cross  sections  for  elastic,  dissociative,  dissociative  with 
attachment,  ionizing,  and  vibrational  excitations  of  N2  and  O2  are  shown  in  Figure  B.2. 
The  data  on  electron  impact  cross  sections  compiled  in  this  figure  were  taken  from  the 
following  sources:  for  elastic  and  vibrational  collisions  of  electons  with  N2  and  O2 
[Phelps,  1985],  for  the  dissociative  attachment  of  electrons  to  O2  [Schulz,  1962],  for  the 
ionization  [Mark,  1975],  for  the  dissociation  from  [Tsang  et  al.,  1991], 

The  specific  value  of  the  rotational  constant  of  N2  is  2.48  x  10-4  eV=  2.9°  K  and  of 
02  is  1.79  x  lO"4  eV=  2.1°  K  [Herzberg,  1950], 

Also  [Murphy,  1988]  is  a  good  general  reference  concerning  electron  collision  cross 
sections  for  O2  and  N2. 
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Electron  Impact  Energy  (eV) 

Fig.  B.l.  Electron  impact  cross  sections  for  optical  emissions.  References  for  each 
particular  line  are  given  in  the  text. 
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Fig.  B.2.  Electron  impact  cross  sections.  References  for  each  particular  line  are  given 
in  the  text. 


c 

Numerical  Algorithms 


C.l  THE  BOLTZMANN  KINETIC  EQUATION 


In  the  numerical  solution  of  equations  (2.4)-(2.6)  we  take  into  account  the  fact  that  the 
total  collision  frequency  u  of  electrons  with  neutrals  is  much  larger  than  the  typical  char¬ 
acteristic  frequency  of  the  EM  pulse  ( d/dt )  and  neglect  the  time  derivative  in  equation 
(2.5).  To  support  this  assumption  we  note  that  even  at  100  km  altitude  the  collision  fre¬ 
quency  for  cold  ambient  electrons  is  ~3x  104  Hz  (see  for  example  Figure  1  in  [ Taranenko 
et  al.,  1992]),  which  is  still  larger  than  10  kHz,  the  value  of  the  central  frequency  of  the 
pulse  used  in  our  calculations.  Also,  the  collision  frequency  increases  substantially  dur¬ 
ing  heating  and  exponentially  increases  with  decreasing  altitude,  so  that  this  assumption 
is  fully  valid  for  the  problem  at  hand  over  the  altitude  range  of  interest. 


After  neglecting  the  time  derivative  in  equation  (2.5),  f\  can  be  expressed  in  a  straight 
forward  way  as  a  function  of  /o  and  such  an  expression  can  be  substituted  into  (2.4) 
to  give  us  an  equation  describing  the  evolution  of  /o  under  the  influence  of  the  electric 
field.  For  example,  for  E  ||  B0  the  procedure  described  above  gives 


dfo  _  e2E2  d 

(v2  dfo\ 

i+  i  d ' 

dt  3 m2v2  dv 

U  dv) 

2v 2  dv  . 

v2Svd  ~  +  vf0)  +  Soi„  (C.  1 ) 


Without  S’om  this  equation  is  a  Fokker-Plank  equation  which  appears  in  many  branches 
of  physical  sciences  [Risken,  1984].  Hence,  one  of  the  ways  to  treat  (C.l)  is  to  consider 
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it  as  a  combination  of  the  Fokker-Plank  equation  and  an  operator  Som  that  transports 
electrons  to  lower  energies  due  to  inelastic  loss  processes  at  each  time  step  (assuming 
discrete  representation  for  the  time  and  velocity  variables).  In  our  computer  model  we 
used  the  Chang  and  Cooper  [1970]  differencing  scheme  to  represent  the  Fokker-Plank 
part  of  equation  (C.l). 

Another  way  to  solve  equation  (C.  1 )  is  to  convert  it  to  the  Boltzmann  transport  equation 
for  electrons  in  energy  space  and  represent  it  in  discrete  form  [Rockwood  and  Green, 
1980]  as 


(C.2) 


where  n*  is  the  number  of  electrons  per  unit  volume  per  unit  energy  in  the  energy 
interval  (ejt+i  —  £*).  The  matrix  elements  C/jt  are  constants  (for  a  given  electric  field) 
which  represent  the  rate  of  transfer  of  electrons  from  energy  £/  to  £*  due  to  elastic  and 
inelastic  collisions. 


It  turns  out  that  both  of  the  two  approaches  are  basically  the  same  since  in  the  discrete 


representation  velocity  derivatives  describe  transport  of  electrons  between  adjacent  grid 
cells. 

A  FORTRAN  computer  code  for  solving  equation  (C.  1 )  can  be  obtained  from  the  VLF 
group  of  STAR  Laboratory  of  Stanford  University. 

C.2  MAXWELL  EQUATIONS 

To  describe  the  propagation  of  EM  fields  we  numerically  solve  Maxwell’s  equations  (2.1). 
Each  of  the  two  Maxwell  equations  can  be  treated  as  a  parabolic  equation  in  discrete  space 
and  time  meshes  and  any  of  many  different  numerical  methods  for  such  equations  can  be 
applied.  A  good  summary  of  such  methods  can  be  found  in  Table  3.1  of  Potter  [1973]. 
In  our  realization  we  used  the  Lax-Wedroff  method.  The  mapping  of  the  electric  field 
amplitudes  to  the  electric  current  density  was  done  by  solving  equations  (2.4)-(2.6)  for 
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the  electron  distribution  function  in  the  way  described  above.  The  multi-slab  technique 
that  uses  the  steady  state  solution  of  the  Boltzmann  equation  is  described  in  Section  2.4. 

A  FORTRAN  computer  code  for  solving  Maxwell’s  equations  (2.1)  as  applicable  to 
this  particular  problem  can  be  obtained  from  VLF  group  of  STAR  Laboratory  of  Stanford 
University. 


D 

Quenching  and  Emission  Rates 


This  appendix  consists  of  a  table  of  quenching  and  emission  rates  given  on  the  next 
page.  Each  reaction  in  the  Table  that  contains  an  hu  on  its  r.h.s.  represent  optical 
emissions  while  others  represent  quenching  reactions.  For  convenience  of  readers  the 
optical  transitions  are  summarized  as  follows:  6300  (6364)  A  {XD  — »  3P  transitions), 
and  the  green  5577  A  (XS  — »  lD  transitions)  lines  of  O;  the  blue  to  UV  2nd  positive  band 
(C3 7ru  — *  B*irg  transitions)  and  red  to  IR  1st  positive  band  (P3^  — ►  v43E„  transitions) 
of  N2;  the  blue  1st  negative  band  (P2£*  — < ►  A'2L^  transition)  and  IR  Meinel  band 
(j427tu  — >  X2l.g  transition)  of  Nj  ions;  and  the  red  1st  negative  band  (64Z~  — »  a47ru)  of 
Oj  ions. 
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APPENDIX  D:  QUENCHING  AND  EMISSION  RATES 


TABLE  Dl.  Quenching  and  Emission  Rates 


Reaction 

Rates  (cm3  s  1  or  s  1 ) 

0 ('DH  O+hv 

9.34  x  10"3* 

0(,D)+N2  o+n2 

2.0  x  10“ 11  exp(107.8/T)* 

0(1iD)+02  — ►  O+O2 

2.9  x  10-11  exp(67.5/T)* 

(X'DHO^  O+O 

8  x  10"12* 

O('S)— ►0(1D)  +  hv 

1.07* 

0 (lS)~*  O+hv 

0.0444* 

0(1S)+02  ->  0+02 

4.9  x  10"12exp(— 885/T)* 

0('S)+0— ►  O+O 

2  x  10" 14* 

N2(C37ru)  ►  N2(B\g)  +  hv 

2xl07* 

N2(C37ru)+02  — ►  o2+n2 

3  x  10-10t 

N2(A3I+)  +  hv 

1.7xl05t 

N2(537r5)+N2  -►  N2+N2 

io-ut 

N|(52i:)^  N+2(X21 +g)  +  hv 

1.4x  1 06$ 

N^(B21+J+N2  -  n2+n2 

4  x  10-10t 

N£(A2ttu)^  N+(X2Z +g)  +  hv 

7xl04t 

n;(b2i;)+n2  -+  n2+n2 

5  x  10-*°$ 

02(b4l,g  )— ►  02(a4nu)  +  hv 

8.5x  104t 

0|(64Z7)+N2  -+  o2+n2 

2  x  10-10$ 

*  Torr  et  al.  [1990] 

$  Vallance  Jones  [1974]  (p.  115,  119) 
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